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Chapter 1 



Introduction 



Ordinary hadronic matter, i.e. matter made of protons and neutrons, is composed of 
elementary particles called quarks and has quantum numbers as dictated by the quark 
model [1]. The quark model is a non-dynamical theory that classifies the different hadrons 
in terms of their quark content. This model allows two families of hadrons: mesons, 
formed by a quark and an antiquark and baryons, formed by three quarks. States found 
experimentally but not describable by the quark model are called exotic states. 

Quantum Chromodynamics (QCD) is the theory that describes the strong interaction 
[2], which is responsible for the formation of hadrons as bound-states of quarks. It is a 
quantum gauge-field theory that has quarks and gluons as the elementary degrees of free- 
dom (the quantization process additionally introduces unphysical auxiliary fields, called 
ghosts). The charge corresponding to the strong interaction in QCD is called color charge. 
It is an yet unproven statement (and so far an extraordinarily well-established experimen- 
tal fact) that all physical states must have no net color charge; this is the most intuitive 
definition of color confinement. 

QCD was built upon the understanding of hadrons via the quark model and was 
intended to provide a dynamical description of the formation of hadrons as states composed 
of quarks and bound by gluons. However, the consensus that QCD provides a correct 
picture of strong interactions is due to its success at describing high-energy processes. The 
reason for this is that at high-energies QCD becomes a weakly-coupled theory and therefore 
perturbative methods can be applied in this regime. The calculation of hadron properties 
developed, instead, at a lower pace. The reason for this is two-fold: at low energies QCD 
becomes a strongly-interacting theory and, on the other hand, bound-state formation is 
an essentially non-perturbative phenomenon. In conclusion, perturbative techniques no 
longer apply in the low-energy regime and new methods had to be developed. 

The challenge of low-energy QCD calculations is to understand how hadrons emerge 
as the physical degrees of freedom, out of the elementary degrees of freedom of the theory, 
namely quarks and gluons (and ghosts). At the phenomenological level, this amounts to 
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calculate hadron properties (mass, radius, shape, etc.) from the QCD Lagrangian 



where q, A and c are the quark, gluon and ghost fields, respectively. Latin indices represent 
color (in the adjoint representation; quarks carry a color index in the fundamental repre- 
sentation, which is ommitted here for simplicity) and A are the Gell-Mann matrices. This 
expression defines the renormalization constants, which are not completely independent 
but have to fulfill the following relations 



Z\f — Z g Z2Z^ 2 , Z\ — Z g Z^ 2 , Z\ — ZgZ^Z^ 1 , Z^ — Z 2 g Z\ (1.2) 



which are consequence of the Slavnov- Taylor identities. 

The spatial distribution of some of the baryon properties, such as mass or electric 
charge, is of especial relevance to understand low-energy QCD dynamics, since they probe 
the details of the the quark-quark and quark-gluon interactions. In particular the ques- 
tion arises whether the shape of baryons deviate from sphericity. For this reason, the 
electromagnetic properties of the nucleon and the Delta resonance have been subject to 
extensive experimental research. The evolution of the nucleon electromagnetic form fac- 
tors with the photon momenta is now very well known experimentally (see, e.g. [3, 4] for 
recent reviews). On the other hand, the experimental information for the Delta electro- 
magnetic form factors is very limited due to the short lifetime of this resonance. Most of 
this information comes indirectly from the study of the jN — > A transition [5-12]. The 
only direct information on the Delta electromagnetic properties is limited to the A ++ and 
A + magnetic moments [13, 14], but with large errors. A more precise measurement of the 
A+ magnetic moment is expected to be performed at MAMI [15]. With this state of af- 
fairs, any theoretical calculation of the Delta electromagnetic form factors which can lead 
to model-independent statements constitute a prediction. We will address this problem in 
Chapter 3. 

Traditionally, hadron properties have been studied by modeling QCD with effective 
degrees of freedom. For instance, constituent quark models (see e.g [16, 17] and refer- 
ences therein) describe hadrons as bound-states of effective quarks, whose mass includes 
dynamical effects (the typical constituent quark mass is one-third the nucleon mass), and 
use a wide variety of, more or less complicated, interaction potentials among the quarks. 
They have been very successful in describing hadron spectra as well as other properties 
such as electromagnetic form factors. The main problem of these approaches is that it is 
not clear whether they capture all features, or any, of QCD dynamics. 

Lattice QCD methods provide non-perturbative calculations of hadron properties using 
the fundamental degrees of freedom of QCD (see e.g. [18] and references therein). In many 
cases they almost have the status of theoretical experiments, when no experimental data 




+ Z 3 c a d 2 c a - Z 1F ig q^-qAl - Z l9 f abc (d,A a v ) A^A* 
+ Z±\g 2 f abe f de A«A b u A^Ai + Z l9 f abc c a d,(A c ^) , 



(1.1) 
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for the situation of interest is available (for example, in Chapter 2 we calculate the mass of 
triple-charm and triple-beauty baryons which have not been observed yet and, therefore, 
we can only compare to lattice predictions for their masses). Lattice approaches have, 
nevertheless, some drawbacks. First of all, most of the calculations are performed in the 
so-called quenched approximation, in which the quantum fluctuations involving quarks are 
neglected. Also, common to all lattice calculations, are the problems of discretization and 
finite-volume effects. Moreover, for technical reasons, the calculations are performed at 
unphysical quark masses (or pion masses) and some procedure to extrapolate the results 
to the physical mass must be defined. With the increase of computing power, however, 
these limitations are rapidly diminishing. An intrinsic problem of lattice methods is that 
it is difficult to unravel how bound-states are formed since, by the very nature of lattice 
calculations, they include all possible quantum correlations among quarks. 

Hadron physics should, in principle, be describable directly from QCD as a contin- 
uum quantum field theory. A rigorous, and systematically improvable, method to extract 
hadron properties from the QCD Lagrangian (and widely used to extrapolate lattice cal- 
culations to the physical pion mass) is chiral perturbation theory (for a pedagogical intro- 
duction see, e.g., [19, 20]). However, this approach is limited to the light quark and low 
momentum region. 

A complete description of a continuum quantum field theory, and in particular of QCD, 
is given when all the (infinitely many) Green's functions of the theory are known. Func- 
tional methods (e.g., Functional Renormalization Group and Dyson-Schwinger equations) 
provide such a description. In particular, Dyson-Schwinger equations (DSEs) [21, 22] 
constitute an infinite set of coupled, non-linear integral equations for the full Green's 
functions of the theory (in a quantum field theory one must distinguish between bare 
Green's functions, which are derived directly from the Lagrangian and are purely classi- 
cal, and the full or dressed Green's functions, which include all quantum effects). These 
are the same Green's functions that one studies in lattice QCD, of which therefore DSEs 
offer a complementary approach. An interesting feature of these equations is that they 
can be solved exactly in the infrared-momentum region, thus providing insight into the 
non-perturbative regime of the theory. However, any feasible solution of DSEs in a general 
momentum range requires to truncate the system to a finite number of equations (for a 
review see, e.g,. [23, 24]). 

Hadrons are encoded in those Green's functions. In general, a bound state in quantum 
field theory corresponds to a pole in the full Green's function describing the evolution of 
the constituent particles (see e.g. [25]). This pole cannot be obtained from a perturbative 
expansion in Feynman diagrams but requires a non-perturbative treatment. In the DSE 
framework, bound-states are described by generalized covariant Bethe-Salpeter equations 
(BSEs). They can be defined from the corresponding Green's function and its defining DSE 
by a Laurent expansion around the bound-state pole. The bound-state is now represented 
by the so-called Bethe-Salpeter amplitudes, which are related to the residue of the Green's 
function at the pole and are the solutions of the BSEs. These equations require, as an 
input, some of the Green's functions obtained from the full DSE system. Naturally, if 
a truncation of the DSEs is performed, a consistent truncation in the BSE of interest is 
required. Nevertheless, once such a truncation scheme is fixed, BSEs provide information 
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both about the particle spectrum and their internal composition. 

Mesons have been thoroughly studied within the DSE/BSE framework, mostly using 
the simplest of the truncations, so-called Rainbow-Ladder truncation. In this truncation 
scheme, of all the possible interactions between the two quarks forming the meson, only a 
single dressed gluon-exchange is taken into account. Moreover, the quark-gluon interaction 
vertex and the gluon propagator are modelled (and restricted to depend only on the mo- 
mentum of the exchanged gluon) and only the quark propagator is solved self-consistently 
from its DSE (see, e.g., [26-28]). There are also some studies about the role of beyond 
Rainbow-Ladder effects [29-31]. 

Baryons, being a three-body system, are much more complicated and therefore have 
been less studied so far but, nevertheless, significant progress using also the Rainbow- 
Ladder truncation has been achieved recently [32-35] . Since the Rainbow-Ladder trunca- 
tion entails the choice of a model for some dressing functions, it is not entirely transparent 
what features of the calculation are due to the modeling and which ones to the truncation 
itself. The goal of this thesis is to make some model-independent statements about the 
study of baryon properties using covariant Bethe-Salpeter equations within the Rainbow- 
Ladder truncation scheme. To do this we perform the calculations using two models very 
different in nature. 

Another interesting aspect of bound-states in QCD is to study the nature of exotic 
hadrons. Since QCD is a non-abelian Yang-Mills theory, the gauge bosons have self- 
interactions. For this reason, QCD predicts that there should exist bound states formed 
by gluons only, so-called glueballs [36]. Although there is an intense experimental effort 
to discover glueballs, for the moment there is no direct evidence of them. In Chapter 4 we 
propose a BSE for a two-gluon system. To derive it, we start from the DSE for the four- 
gluon Green's function, which describes the propagation of two gluons in spacetime, and 
select the diagrams which would develop a pole if a bound-state is formed. A consistent 
solution of this equation, which is beyond the scope of this thesis, would require the 
knowledge of several full, i.e. dressed, Green's functions, namely the gluon propagator, 
and full three-gluon and four-gluon vertices. These Green's functions can, in principle, be 
obtained by solving the corresponding DSEs. 

In the next section we give a brief description of how to derive the Dyson-Schwinger 
equations from a given Lagrangian 

1.1 Dyson-Schwinger and Bethe-Salpeter equations in a nut- 
shell 

As already mentioned above, a bound state in quantum field theory corresponds to a pole 
in the corresponding Green's function. Those poles cannot appear in a perturbative series, 
but are an essentially non-perturbative phenomenon. The residue of the Green's function 
at the the pole allows to define the Bethe-Salpeter amplitudes (see e.g. [25]). 

The steps to write down a relativistic equation for bound states can be summarized 

as: 

• Derive a non-perturbative equation describing the relevant Green's function. 
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• Examine the structure of the equation to find those terms which will develop a pole 
when a bound-state is formed (pole ansatz). 

• Performing a Laurent expansion of the Green's function around the pole one finds 
an homogeneous equation for the Bethe-Salpeter amplitudes, the Bethe-Salpeter 
equation, and a normalization condition. 

In the rest of this section we shortly describe the first point. The other two will be 
developed, when necessary, along the thesis. 



1.1.1 Equation for the Green's function: Diagrammatic derivation 

In this subsection we describe a diagrammatic derivation of an equation for a given Green's 
function (for details see e.g. [37]). Equations of this type are commonly referred to as 
Dyson equations. Although the formal derivation given in the next subsection is more 
powerful and does not rely on any kind of perturbative expansion, the steps described 
here can provide a more intuitive, or at least complementary, picture of how the different 
interaction terms appear. 

A 2n-points Green's function is given by the (physical-) vacuum expectation value of 
the time-ordered product of 2n Heisenberg field operators 



G 



a 1 ...a n ;a 1 . 



a , n {x X ...X n ; X [... x'J = (0 out \T [Vax^l) • • • K) • • • ^a'S X n)\ , 

(1.3) 



where the generic indices a represent all possible indices carried by the fields. By using 
the iterative expansion of the time evolution operator (see, e.g [38]), the above expression 
can be rewritten in terms of interaction-picture operators: 



Gai...a n ;a' 1 ...a' n (^l ■ ■ ■ x n] x \ ■ ■ ■ x n) — 

d 4 y 1 ...d 4 y k (0\T[^ :ai (x 1 )...^ an (x n ) 

'4), a 'S x i)--^la'y n )ni{yi)..M I {y k )] |0> , (1.4) 



where the k = case represent the free propagation of fields. This perturbative expression 
is suitable to study, for instance, scattering processes. However, if the system develops 
bound states the Green's function will have poles. A pole in the Green's function will 
never appear by summing any finite number of diagrams in this expansion, but instead 
the whole series, or an infinite subset of it must be considered. 

Applying Wick's theorem to the elements on the right-hand side of (1.4) one generates 
all possible diagrams (see Figure 1.1). The (infinite) sum of all ^-particle-irreducible 
connected terms (with 1 < £ < n) is called irreducible ^-particle interaction kernel —iK^ 
and the sum of all one-particle terms gives the full propagators S. All reducible diagrams 
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G <3> - = + 



+ 



+ 



+ 



Figure 1.1: Some examples of diagrams that would be generated by Wick 
contractions in Equation (1.4) for the case of a 6-points Green's function. The 
full series contains all possible one-particle- (e.g. second diagram), two-particle- 
(e.g. third and fourth diagrams) and three-particle-irreducible (e.g fifth dia- 
gram) terms. 



can be generated by iteration of the irreducible ones, using the full propagators as internal 
lines. Defining the modified kernels as 

K®= Yl K^ S-\ v S-\ , (1.5) 

n—l 

cyclic perm. 
£ elements 

we can introduce a single interaction kernel 

K= ^ (1.6) 

1=2. ..n 

and write an inhomogenous integral equation for the Green's function 

Ga 1 ...a n ;a' 1 ...a' n i x l ■ • ■ x n'i x \ ■ ■ ■ x n) = ^aia[ ( x li x l) ■ • • ^a„a' n ( x m x n ) 

-i / d 4 yi...d A y n S aibl (x 1 ,y 1 )...S anbn (x n ,y n ) 



x J d 4 y[ . ..dy n K bl _ hn . Vi ... Vn (y 1 . ..y n ;y[ . . . y^G^...^^..^^! • ■■Vn^'i ■ ■ - x 'n) ( L7 ) 
or, symbolically 

G = G - iG KG (1.8) 

where Go represents the product of full propagators. This is the Dyson equation for the 
Green's function G. 

It is convenient to work in terms of the scattering matrix T, defined by amputating all 
incoming and outgoing legs 

G = G + G TG (1.9) 

and substituting in (1.8), one gets 

T = -iK-iKG T. (1.10) 

The advantage of having an integral equation for the Green's function is that, even making 
approximations for the interaction kernel K, one is considering an infinite set of diagrams 
and therefore is useful to study bound states. 
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1.1.2 Dyson-Schwinger equations: Formal derivation 

In this section we describe the derivation of Dyson-Schwinger equations from the action 
of a given quantum field theory. These equations relate a particular Green's function 
with higher order Green's functions which, in turn, fulfill their own Dyson-Schwinger 
equations. They constitute an infinite and coupled system of integral equations for the 
Green's functions of the theory which do not rely on any perturbative expansion. 

A quantum field theory is defined (in Euclidean spacetime) by the generating functional 

z\j\ = j v $ e -sm+s*xux)* a {x) = Gii inJii . . . Jin j (i.n) 

where & a are generic fields, J a are the sources associated to these fields and we defined 
here the full Green's functions G^...^ as the moments in an expansion in terms of the 
sources, with the indices ij denoting any possible discrete index of the field as well as 
spacetime variables. 

Assuming that the functional integration J D$ is defined such that the integral of a 
total derivative vanishes, we obtain 

*- m+J '*' ~ I " ^ (I 3 -*) = (I s " *) w ■ 

(1-12) 

where the brackets in the last term denote the vaccuum expectation value with non-zero 
sources. We can take the term in parentheses out of the integral by noting that, acting on 
Z[J], we can make the identification <£j — > 5/SJi. We have then the following identity 



8S 



+ Jt Z[J] = . (1.13) 



This equation provides the seed to generate all Dyson-Schwinger equations for the full 
Green's functions, upon taking appropriate functional derivatives with respect to the 
sources and afterwards setting them to zero. 

As an example, we will derive the Dyson-Schwinger equation for the quark propagator 

S(x-y) = (0\q(x)q(y)\0) , (1.14) 

since it is the essential element in any covariant bound-state calculation for mesons and 
baryons. We start with a particular case of (1.12) 

/<5 ^-^)\ =0, (1.15) 



\ Sq(a 



[J] 



with q(x) the anti-quark field and j(x) its source. Taking a further derivative with respect 
to j(y) we obtain 

(tfH,r 5,4,(1 -*>• (116) 
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Using the explicit expression of the Lagrangian (1.1) to perform the derivative 5Sqcd/$Q 
and afterwards setting the sources to zero, we obtain the quark-propagator DSE 

8W(x -y) = Z 2 {-$> + Z m m)S(x - y) 

- igZ lf J d 4 zd 4 z'5^(x ~ z)S^(x - z') 7 ^( q (z)q(y)A«(z')) , (1.17) 

where the full quark-gluon Green's function (q(z)q(y)A a ^(z')) can be written as a proper 
(one-particle irreducible) quark-gluon vertex F^(x,y,z) with two quark propagators and 
one gluon propagator attached. In momentum space, the quark DSE reads 

/d 4 k t a 
- lfl S(k)r a u (k,p)D IIU (q) . (1.18) 



Chapter 2 

Baryon Masses 



In this chapter we calculate the masses of spin-!/2 and spin- 3 /2 ground-state baryons using 
a covariant three-body Bethe-Salpeter equation, also called covariant Faddeev equation. 
This equation assumes that three-body-irreducible interactions are negligible. Moreover, 
of all possible two-body-irreducible interactions only a single dressed-gluon exchange is 
kept. This truncation induces the necessity of a model for the quark-gluon interaction. 
Using two models, very different in nature, for this interaction we intend to make model- 
independent statements. We will find that this simple setup gives a good description of 
baryon masses, up to a 10% accuracy. 



2.1 Covariant three-body Bethe-Salpeter equation 

The evolution of a three-quark system is encoded in the six-point Green's function G^ 
or equivalently in the six-point scattering matrix . As explained in previous chapter, 
both are described by the corresponding Dyson equations (see Figure 2.1) 

G(3) = G (3) _ iG f) K G (3) 

T (3) = -iK-iKG^TW ( ' j 

where, using the definitions in (1.5), the interaction kernel K can be decomposed into 
four terms and (with a = 1, 2, 3 denoting the spectator quark), containing only 
three- and two-particle irreducible graphs, respectively. 

We already discussed in previous chapter that when the three quarks form a baryon, 
the scattering matrix (and the Green's function) develops a pole at P 2 = —M 2 , where P 
is the total momentum of the bound-state 

P = Pl+P2+P3, (2.2) 

with pi the quark momenta, and M the baryon mass. On the baryon mass-shell we write 
the scattering matrix as 
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G (3) 






-iK 






G (3) 












> 








iK 



J(3) 



Figure 2.1: Dyson equation (2.1) for a three-particle system. The system 
is equivalently described by the Green's function or the scattering ma- 
trix T^ 3 \ Lines with blobs represent fully dressed propagators and K is the 
interaction kernel. 



which defines the Bethe-Salpeter amplitude and its conjugate ^f. The factor C s , which 
depends on the baryon spin s, is 

s = \ : C s = 2MA+(P) , (2.4) 



s = I : C s = 2MP^(P) , (2.5) 



where we introduced the positive-energy Dirac projector A + and the Rarita-Schwinger 
positive-energy projector P^ u 

k + (P)= l -{t + f) , (2.6) 

Pnp) = A + (P)(V-i 7 £ 7 ^ , (2-7) 

with 7^ = Tp U, y u , Tp U the transverse projector 

T^{q) = V - ( 2 - 8 ) 

and the hat denotes a unit vector 

t>" = -= . (2.9) 

Inserting this expression in (2.1) and neglecting regular terms, we can write an equation 
for "if, the three-body Bethe-Salpeter equation (see Figure 2.2) 

3 

= -iK^ G ( q ] y + Y, ~ iK (a) G o ] * » ( 2 - 10 ) 

o=l 

where we used (1.5) to introduce de modified kernels K. 

The nature of the BS amplitudes \& depends on the baryon of interest (see Appendix 
A. 2): for spin-i/2 baryons it is a rank-4 Dirac tensor and for spin- 3 /2 baryons it is a mixed 
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tensor with four Dirac and one Lorentz indices (and on top of this they have flavor and color 
structure) . Its general structure can be determined imposing only Poincare covariance and 
parity invariance. To solve this equation, one needs to specify the interaction kernels and 
the full quark propagator. 




Figure 2.2: Pictorial representation of the three-body Bethe-Salpeter equa- 
tion (2.10). Lines with blobs represent fully dressed propagators. The kernels 
and K ^ contain only three- and two-particle irreducible graphs, respec- 
tively. 



2.2 Quark propagator and Rainbow truncation 

The full quark propagator is an essential element for any covariant bound-state calculation 
in QCD. The general structure of this propagator is 

(S- 1 )i B (p)=A(p 2 ) (i^ + M(p 2 )) ab 8 AB (2.11) 

where 1/A(p 2 ) is the quark wave-function renormalization and M(p 2 ) is the quark mass 
function. The indices a, b and A, B are Dirac and color indices, respectively. 

In this work, the quark propagator is obtained by solving the quark Dyson- Schwinger 
equation (1.18), 

(S- 1 )i B (p) = Z 2 (ip' + Z m m q ) ab 




where q = k — p is the momentum of the exchanged gluon and here r, s represent color 
indices in the adjoint representation. The color traces can be worked out (see Appendix 
A. 3) to give (N 2 — 1) /(2N C ), with N c the number of colors. The renormalization constants 
Zif(fi 2 , A 2 eg ), ^2 (a* 2 , A 2 e9 ) and Z m (fj?, A 2 eg ) depend on the renormalization scale p, and 
on an ultraviolet scale A reg , required to regularize the integral. The renormalized mass 
m q is related to the bare mass mo by mo = Z m m q . Equation (2.12) is solved imposing 
the renormalization conditions 



A(p 2 = v 2 ) = 1 , 
M(p 2 = n 2 ) = m , 



(2.13) 
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q 




k=p-q 



Figure 2.3: Pictorial representation of the quark Dyson-Schwinger equation 
(2.12). Blobs represent fully dressed propagators or vertices. 

at a sufficiently large [i and for a fixed value for mo- 

A diagrammatic representation of this equation is given in Figure 2.3. In Figure 2.5 
we illustrate the dynamical generation of quark mass when the quark DSE is solved using 
appropriate models for the quark-gluon vertex, as explained in the next sections. We also 
show how this feature disappears if the interaction is not strong enough. 

The Green's functions required to solve the quark propagator DSE (2.12) are the full 
gluon propagator D fiu (q) and the full quark-gluon vertex r^(p, q) (we omit here all other 
indices). Slavnov- Taylor identities restrict the full gluon propagator, in Landau gauge, to 
be 

D-M = S rs (<W - qj f) ^f 1 (2-14) 

with only one unknown dressing function Z{q 2 ). On the other hand, using symmetry 
arguments, one can decompose the full quark-gluon vertex into 12 linearly independent 
Lorentz covariants T-~ 

12 

r"(p, q) = Y, <?>P ■ 2) T i>> ?) ( 2 - 15 ) 

i=l 

where fi are Lorentz-invariant dressing functions. Eight of these twelve components are 
purely transverse. 

A completely consistent solution (that is, using a gluon propagator and quark-gluon 
vertex obtained as solutions of their own DSEs) is not possible in a general momentum 
range* and some truncation scheme must be chosen. The simplest possibility in the context 
of the quark DSE is to use for the full quark-gluon vertex the tree-level vertex 7^ times a 
dressing function. More specifically, if we single out the tree-level part in the full quark- 
gluon vertex as 

Y> t {p,q) = Z xr f + M t , (2.16) 

with all other structures gathered in A' 4 , one can truncate the full vertex by projecting it 
onto the tree- level part and limiting it to depend only on the gluon momentum q, 

T^p,q)^(Z lf + A(q 2 ))^ . (2.17) 

This is known as the Rainbow truncation of the quark DSE. The flavor and color parts 
of the vertex are the tree- level ones. Within this truncation, only the dressing functions 



*It is possible, nevertheless, to solve the full set of DSEs in the deep infrared momentum-regime, but 
this regime is not expected to be of relevance to hadron physics. 
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Z(q 2 ) and A(g 2 ) have to be fixed. It has been customary in previous phenomenological 
studies of hadron properties in the DSE/BSE approach to combine all scalar dressings 
into a renormalization-group invariant ^ effective running coupling, a e ff(q 2 ) 

Z if D £W(<?) r -(^) -»• Z l T AQ) ^ ■ ( 2 - 18 ) 

If, instead, one wishes to draw a distinction between the gluon and the quark-gluon vertex 
dressings, it can be useful to model them separately. This distinction can be important 
because the gluon propagator is by now well-known from both lattice studies and other 
functional approaches. 

A truncation of the full quark-gluon vertex means that, in a diagrammatic expansion 
of the vertex (or skeleton expansion), only one (or several) of those terms is taken into 
account. Therefore, once such a truncation for the quark DSE is chosen, a consistent 
truncation of the interaction kernel in the BSE must be defined, as described in the next 
section. 

As will be explained later in this chapter, the momenta of the internal quark prop- 
agators in a BSE are complex. The numerical techniques used to solve the quark DSE 
in the complex plane have been described, for instance, in [39, 40]. In this work we use 
the method of [40] and we refer the reader to this article for details. It is worth noting 
that the resulting quark propagator shows complex conjugate poles in the complex plane. 
Whether this is an artifact of the truncation or not is not clear since this feature is found 
in more complex truncation schemes [40, 41]. In any case, the appearance of these poles 
results in a limitation of the maximum bound-state mass one can study (see [39]). 



2.3 Rainbow-Ladder truncation 

The QCD Lagrangian features an approximate SUjt(Nf) £g> SUi(Nf) symmetry, known 
as chiral symmetry (which becomes exact when quarks are considered massless). This 
symmetry would imply, for example, that the n and a mesons or the p and a\ mesons are 
degenerate in mass, which is not the case. In nature, this symmetry is spontaneously bro- 
ken due to the dynamical generation of quark mass and the associated Nambu-Goldstone 
bosons are the pions. 

If mesons are studied using a covariant quark-antiquark BSE, and for the quark prop- 
agators one takes the solutions of the quark DSE, it has been shown [42, 43] that the 
crucial relation that ensures a correct implementation of chiral symmetry and its dynam- 
ical breaking is the axial-vector Ward-Takahashi identity (AxVWTI) 

-i P,T^(p-P) = S-\p + + 1 J-^S-\p + P/2) , (2.19) 

^This can be seen taking into account that, under a change of the renormalization scale, g scales as 
the quark-gluon vertex as and the gluon propagator as I/Z3, and using the Slavnov- Taylor 

identity Z lf = Z g Z 2 Zl /2 
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which relates the axial- vector vertex T^^(p;P) (and, indirectly from its BSE, the quark- 
antiquark kernel) to the quark self-energy. In this equation are the Gell-Mann flavor 
matrices and P and p represent the total and relative quark momentum, 

P = Pi + P2 , 

Pi - P2 (2.20) 

with pi the quark momenta. It can be proven that, in the chiral limit, if (2.19) is fulfilled, 
the dynamical generation of quark mass is accompanied by the appearance of massless 
pseudoscalar bound-states (the pions). 

From (2.19), and given a truncation of the quark DSE, one can find the corresponding 
quark-antiquark interaction kernel. In the case of the Rainbow truncation it corresponds 
to a single dressed-gluon exchange between the quark and the antiquark, with the gluon 
interacting with quarks via a vector coupling. In terms of a e //, it reads 

-iK*~« = 4vr Z\ T^ u (q) i-f i-f (2.21) 

where, as before, the flavor and color parts of the vertex are the tree-level ones. This, to- 
gether with (2.18), constitute what is known as the Rainbow-Ladder (RL) truncation 
of the quark-DSE/meson-BSE system. 

The Rainbow-Ladder truncation provides a flavor-blind and quark-mass independent 
interactions. Corrections beyond Rainbow-Ladder, which may include corrections to the 
effective coupling and the inclusion of additional structures beyond a vector-vector inter- 
action, are expected to change this. An open question is how important are these effects. 
We will see later that a quark-mass dependent interaction is essential to achieve a precise 
description of baryon spectra for all current-quark masses and to observe non-analiticities 
in the current-quark mass evolution of baryon masses due to the opening of decay channels. 

As a side remark, a systematic way to relate truncations in the quark DSE with 
truncations of the BSE kernel such that they preserve (2.19) has been described in [44]. In 
this paper both the quark DSE and the meson BSE are derived from a chirally-symmetric 
2PI effective action so that they automatically fulfill (2.19). Thus, if a chiral-symmetry 
preserving approximation is performed for the effective action, the truncated DSE and BSE 
derived from this action will still fulfill the AxVWTI. As a rule of thumb, if the quark- 
gluon vertex in the quark DSE is truncated in a certain way (that is, keeping only some 
diagrams in its skeleton expansion), a symmetry-preserving BSE kernel can be obtained 
by cutting one internal quark line of those diagrams in all possible ways. 

The goal of this work, however, is the study of baryons within the DSE/BSE framework. 
To this end one needs to fix the three-quark and two-quark irreducible interaction kernels. 
In particular, the two-quark (or quark-quark) kernel is not restricted by (2.19) and one is, 
in principle, free to choose any other truncation scheme. There is, nevertheless, extensive 
literature on the calculation of meson properties using the covariant DSE/BSE approach 
within Rainbow-Ladder (see e.g. [26-28] and references therein). With the idea of having 
a common approach to hadron properties, the Rainbow-Ladder truncation has been also 
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used in baryon studies [32, 33, 45-55]. In this thesis we continue this trend and use (2.21) 
as the quark-quark interaction kernel. 

On the other hand, the study of baryons with covariant bound-state equations was 
performed, until recently, by reducing the three-body problem to a two-body problem 
using the so-called diquark ansatz [45, 46, 56]. It assumes that two-body correlations are 
dominant in baryons and that two of the three valence quarks are bound into an object 
called diquark. This diquark is in turn bound to the remaining quark to form the baryon. 
This approach has been very successful describing baryon phenomenology. 

Moreover, recent calculations of heavy-baryon masses using perturbative non-relativistic 
QCD (pNRQCD) [57] suggest that irreducible three-body interactions contribute only 
~ 25 MeV to the total mass. These results, therefore, give support to the approximation 
of neglecting three-particle irreducible correlations (Faddeev approximation). The 
corresponding three-body Bethe-Salpeter equation is known as covariant Faddeev equation 
(see Figure 2.6), for historical reasons, and in what follows we will refer to the amplitudes 
^ as Faddeev amplitudes. 

2.4 Effective interactions 

As explained in previous section, to completely specify the covariant Faddeev equation one 
needs to model the gluon propagator dressing function and the quark-gluon interaction. 
In fact, this is the only model input of the approach. Therefore, to assess the model- 
independent features of the Rainbow-Ladder truncation, in this work we use two different 
models for the effective interactions. 

The common feature between QCD and both effective interactions is that they repro- 
duce the one-loop behavior of the QCD running coupling at high momentum 

^ - ^ife (2 - 22) 

with 7 m = 12/(lliVc — 2Nf) the anomalous dimension of the quark propagator. The 
models differ in their infrared behavior. 

The first model we use is known as Maris- Tandy model [58, 59] and it has dominated 
hadron studies within Rainbow-Ladder. This dominance is well-earned since this ansatz 
performs very well when it comes to the purely phenomenological calculation of ground- 
state meson and baryon properties. However, this model has no clear connection to QCD 
in the infrared and is, therefore, not entirely satisfactory to gain understanding of the for- 
mation of hadronic bound-states in QCD. On the other hand, with the rapid improvement 
in our knowledge of QCD Green's functions from both lattice and functional approaches, it 
is possible to define different effective interactions which, presumably, capture more faith- 
fully some of QCD's features. Based on this, an effective interaction has been proposed 
in [60]. 

Note that the fact that an effective interaction captures more features of QCD does 
not necessarily mean that it will perform better phenomenologically. This is because 
the interaction is used within a given truncation scheme and, therefore, if one wants to 
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reproduce hadron properties the model has to be tuned to account for the effect of the 
missing contributions. In particular, it has been shown [41] that dynamical quark-mass 
generation is accompanied by the appearance of scalar components in the quark-gluon 
vertex. These components, of course, are missing in the Rainbow-Ladder truncation and, 
therefore, the effective interaction must somehow mimic its effects. This can already be 
seen in Figure 2.5, where dynamical chiral symmetry breaking is missing if the effective 
interaction is too weak. 

In this respect, both models described below are designed to reproduce correctly dy- 
namical chiral-symmetry breaking as well as pion properties at the physical u/d mass. 
This means that they capture beyond Rainbow-Ladder effects at this quark mass. As a 
consequence, both interactions have similar strength at the intermediate momentum region 
~ 0.5 — 1 GeV (see Figure 2.4). To analyze whether the effects beyond Rainbow-Ladder 
are analogous for mesons and for light baryons or, more precisely, whether baryon spectra 
is also well reproduced at this, or at higher, quark mass is one of the goals of this chapter. 




q ! (GeV ! ) q 2 (GeV ') 



Figure 2.4: Left panel: Maris- Tandy effective interaction for different values 
of r\. Right panel: Comparison of the AFW and the MT (for rj = 1.8 models 
using a log-log scale to stress the qualitatively different behavior in the deep 
infrared. 



2.4.1 Maris-Tandy model 

In the Maris-Tandy (MT) model [58, 59] the effective running coupling is given by 

2 \ 2 ,, 2 2ii-'y_ h - P ~i 2 / A t 



27r 7TO (l-e-9V A *) 



« e //(, 2 ) = w )e** + ln[e2 -:r; v a+g2/A ^ )2 ] ■ ( 2 - 23 ) 

which, behaves as (2.22) in the UV and features a Gaussian distribution in the infrared 
(see Figure 2.4) that provides dynamical chiral symmetry breaking. The scale Aj = 1 GeV 
is introduced for technical reasons and has no impact on the results. Therefore, the inter- 
action strength is characterized by an energy scale A, fixed to A = 0.74 GeV to reproduce 
correctly the pion decay constant from the RL-truncated meson-BSE. The dimensionless 
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parameter r/ controls the width of the interaction (see Figure 2.4). Many ground-state 
hadron observables have been found to be almost insensitive to the value of rj around 
rj = 1.8 [28, 33, 51]. This has been used as an argument in favor of the model indepen- 
dence of Rainbow-Ladder results. Instead of pursuing this line of research, we prefer to 
introduce a new, not-related model to evaluate the validity of those assertions. For the 
anomalous dimension we use 7 m = 12/(llNc — 2Nf) = 12/25, corresponding to Nf = 4 
flavors and N c = 3 colors. For the QCD scale Aqcd = 0.234 GeV. 

Note that in the numerical resolution of the quark DSE we employ the Pauli-Villars 
regularization method of the integrals, with a mass scale of 200 GeV. Moreover, for this 
model, we fit the quark masses, at the renormalization scale fi = 19 GeV, to be 3.7, 85.2, 
869 and 3750 MeV for the u/d, s, c, and b quarks, respectively. 



2.4.2 Alkofer-Fischer- Williams model 

The Alkofer-Fischer- Williams (AFW) model [60] is motivated by the desire to account 
for the [Tritanomaly by the Kogut-Susskind mechanism [61]. The effective coupling 
is constructed as the product of the gluon dressing [62, 63] and a model for the non- 
perturbative behavior of the quark-gluon vertex [41], 

t 2\ n( x W y V K ~ 1/2 / ' a + a uv x\~ 10 ( a uv x\~ 25 ° 

The four terms in parentheses are: the IR scaling of the gluon propagator; IR scaling of 
the quark-gluon vertex; logarithmic running of the gluon propagator; and the logarithmic 
running of the quark-gluon vertex. Additionally, the last two are constructed to interpolate 
between the IR and UV behavior. The remaining terms are defined as follows: 

Ac XrV ( 1 1 \ X =q2 / A YMi 

*=Tjr + ..( 1Y» > ^uv = ^ m ( , y = q 2 /A} R , (2.25) 

1+y l + {y-iy \\nz z-lj _ * , A 

z — 9 / A MOM > 

and ao = 8.915/iVc- Here, Aym = 0.71 GeV is the dynamically generated Yang-Mills 
scale, while Amom = 0.5 GeV corresponds to the one-loop perturbative running. The 
IR scaling exponent is k = 0.595353, and the one-loop anomalous dimensions are related 
via 1 + 70 = -2Sq = § N C "fm, with j m = 12/ (11 Nc - 2Nf). We choose Nf = 5 active 
quark flavors at the renormalization point \x = 19 GeV. The constant C = 0.968 is chosen 
such that a e ff runs appropriately in the UV. Finally, Ajr = 0.42 GeV, As = 6.25, and 
Xb = 21.83 determine the IR properties of the quark-gluon vertex and are fitted such 
that the properties of it, K and p mesons are all reasonably well reproduced. The quark 
masses at fx = 19 GeV are 2.76, 55.3, 688 and 3410 MeV for the u/d, s, c, and b quarks, 
respectively. 
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Figure 2.5: Left panel: Mass function M(p 2 ) for chiral and u/d-quark bare 
masses using the MT and AFW models (the precise value of the bare quark- 
mass depends on the model). Right panel: To illustrate the necessity of a 
strong enough effective interaction to account for dynamical chiral symmetry 
breaking, we plot M(p 2 = 0) as a function of the interaction strength A in the 
MT model. The arrow indicates the value of A in our bound-state calculations. 



2.5 Baryon masses from the covariant Faddeev equation 

The three-body Bethe-Salpeter equation (2.10) in the Rainbow-Ladder truncation (see 
Figure 2.6) reads 



+K aa >pp>(k) 5 , a / Q //(fci)5 /3 / /3 //(fc 2 ) ^a"/3"-yx{p {3 \q i3 \P) 



(2.26) 



where we have absorbed the — % factor into the definition K, so that it is defined as in 
(2.21), 

K aa ,^(k) = -AirC Z\ aeff k f ] T, v (k) 7 1, , (2.27) 

where C = —2/3 stems from the traces of the color matrices (see Appendix A. 3). The 
Faddeev amplitudes depend on the quark momenta p\, p2 and ps, but this dependence 
can be reexpressed in terms of the total momentum P and two relative momenta p and q: 



'l-C)P3-C(Pi+i>2), Pi 
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p 



P2 - Pi 
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P3 



P , i-Cp 
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p + CP , 



(2.28) 
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with £ is a free momentum partitioning parameter. The total momentum is constrained 
by P 2 = —M 2 , with M the baryon mass. 

Our calculations are performed in Euclidean spacetime. Using the conventions of 
Appendix A.l it is easy to realize that the arguments of the quark dressing functions, pj, 
are complex. For instance, 

p 2 = (p + £P) • (p + (P) = p 2 - Q 2 M 2 + 2C i Myftf zx . (2.29) 

where z\ = p-P (see Appendix A.l). Throughout this thesis we use the value £ = 1/3. This 
choice allows for a tremendous simplification of the bound-state equation's solution method 
(see Appendix A. 4). Moreover, as mentioned before, the non-trivial analytic structure of 
the quark propagator in the complex plane limits the maximum bound-state available in 
the approach. The choice £ = l /s for the momentum partitioning parameter maximizes 
the accessible bound-state mass range before hitting the poles of the quark-propagator 
[39]. 

The internal quark propagators S depend on the internal quark momenta ki = p% — k 
and ki = pi + k, with k the gluon momentum. The internal relative momenta, for each of 
the three terms in the Faddeev equation, are 




Figure 2.6: Diagrammatic representation of the covariant Faddeev equation 
in the Rainbow-Ladder truncation (2.26). 

The tensor structure of the Faddeev amplitudes VP is described in Appendix A. 2. The 
numerical resolution of the Faddeev equation (for details on this, see Appendix A. 4) is 
simplified if one expands the spin part of the amplitudes in an orthonormal basis {t^} 

V a ^i(p,q,P) = f (i) (p 2 ,q 2 ,z , Zl ,z 2 ) T ^(p,q,P)® FLAVOR® COLOR , (2.31) 

with zq = pr • 5t, z\ = p ■ P and z 2 = q- P and where the subscript T denotes a transverse 
projection with respect to the total momentum P and the hat, as before, denotes unit 
four-vectors. The color and spin parts of the basis are those of the quark model and, 
since the kernel is flavor and color independent, they factor out. We then solve for the 
Lorentz- invariant coefficients f^(p 2 ,q 2 ,zo,zi,Z2). 

From the solution of the Faddeev equation we obtain the bound-state mass and am- 
plitudes. The latter provide information about the internal structure of the baryon that 
can be used to calculate baryon form factors, as shown in the next chapter. 
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2.5.1 Discussion of the results 

In Table 2.1 we summarize our results for the vector-meson, spin-i/2 and spin-3/2 baryons 
at the u/d, s, c and b quark masses, and compare with lattice or experimental results, 
if available. Note that these results differ slightly from those presented in [35] since, in 
connection with the calculation of form factors (see next chapter), we needed to repeat 
the calculations at much higher precision (see Appendix A. 4). 
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Table 2.1: Computed meson and baryon masses (in GeV) for both MT and 
AFW interactions, compared to experiment. Quantities fitted to their experi- 
mental values are indicated by a f . Since the heavy-Omega baryons have not 
been observed yet, we compare to lattice calculations [64, 65] and a recent 
study from pNRQCD [57]. 



Let us analyze first the results at the n/ci-quark mass. As explained before, both effec- 
tive interactions are conceived to reproduce pseudoscalar ground-state meson properties 
at this quark mass and to this end they must include corrections to the quark-antiquark 
kernel beyond Rainbow-Ladder. The calculated mass for nucleon and Delta agrees very 
well with the physical mass; the difference between the calculated and the physical mass, 
and between the result for the two models, is smaller than 1%. Assuming that beyond 
Rainbow-Ladder corrections are the same for both the quark-antiquark and the quark- 
quark kernels, this seems to indicate that irreducible three-body contributions are neg- 
ligible for light quarks (in agreement with the success of the quark-diquark approach to 
baryon properties). The mass of the p-meson is also well reproduced, although the inter- 
action appears slightly too attractive and the calculated masses are about 4 — 8% smaller 
than the physical mass, depending on the model. 

Figure 2.7 shows the evolution of the bound-state masses with respect to the squared 
pion mass. This, by virtue of the Gell-Mann-Oakes-Renner relation [75] is equivalent to 
plotting them with respect the current quark mass; the plot shows results from u/d to s 
quark masses. Our results are compared with lattice data [66-74] at intermediate quark 
masses and with experimental values for the $-meson and f2-baryon at the s-quark mass. 
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Here one sees that the results for different interactions start to differ as the quark mass 
increases. AFW provides always a weaker binding and, for baryons, MT results appears 
to be systematically below the lattice data or the physical Omega mass. Nevertheless, one 
can say that both interactions compare reasonably well with lattice/experimental data 
within a 10% accuracy. Thus, from these results we can establish the accuracy of RL 
calculations of hadron masses to be ~ 10%. Another way of expressing this is that hadron 
masses seems to be determined at ~ 90% by a single dressed-gluon exchange. 




0.1 0.2 0.3 0.4 0.5 
mi (GeV 2 ) 



Figure 2.7: Evolution of the vector-meson and spin-!/2 and spin- 3 /2 baryon 
masses as a function of the pion mass, for the MT and the AFW models. 
Results are compared to lattice data [66-74]. Experimental values are denoted 
by stars. 

However, there is a missing feature in the quark-mass evolution of the A and p masses. 
In our approach, all particles are considered stable. However, A and p can decay into a 
nucleon or into pions, respectively, via strong interactions and therefore have a significant 
width. Such a decay would manifest itself in a non-analytical behavior of the A/p-masses 
as a function of the pion mass when the decay channel would open [52, 76, 77]. Our 
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truncation scheme does not provide such a mechanism. The possibility of a decay can be 
classified as a beyond Rainbow-Ladder effect and, in particular, would lead to a quark-mass 
dependent interaction. 

To get a glimpse of the behavior of the bound-state spectra at even higher quark 
masses, we calculate the J/tp and T vector- mesons, and triple-charm and triple-beauty 
Omega- like baryon masses. The results are shown in Table 2.1. For the Omega-like baryons 
the are no experimental values and we compare with lattice calculations [64, 65] and a 
perturbative non-relativistic QCD calculation [57] . At the charm quark-mass the trend of 
both interactions is maintained; namely, vector-meson masses are overestimated whereas 
for the Omega-like state the MT and AFW models underestimate and overestimate the 
lattice result, respectively. Also, we can say that both models agree between each other 
and with lattice/physical values within ~ 10%. The result at the beauty quark-mass 
is surprising, since both models give essentially the same result. The reason for this is 
unclear at the moment and deserves further study. 



Chapter 3 



Spin-3/2 electromagnetic form 
factors 

In the previous chapter we calculated baryon and vector-meson masses using two different 
effective interactions for the Rainbow-Ladder-truncated three-body Bethe-Salpeter equa- 
tion. We quantified the model-independent accuracy of Rainbow-Ladder results to be 
~ 10%. We extend here this study to the electromagnetic properties of spin-3/2 baryons 
with spacelike photon momentum. On one hand we try to confirm the validity of our 
estimate for the accuracy of Rainbow-Ladder calculations. On the other hand, due to the 
absence of experimental data and the limitations of theoretical calculations on Delta and 
Omega electromagnetic form factors, our results constitute qualitative predictions of the 
shape of these baryons. 

3.1 Coupling photons to baryons 

A procedure to couple a gauge field to a Green's function described by integral equations, 
in a way such that gauge symmetry is preserved and there is no overcounting of diagrams 
has been described in [78, 79] and named by the authors as gauging of equations. Here 
we reproduce the main steps to derive an expression for a baryon coupled to an external 
electromagnetic field. 

The idea is that, starting with a 2n-points Green's function 

G^\x x ...x n ;x[... x' n ) = (0 out \T [^(xi) . . . V n (x„)^ Vi) • • • \0 in ) , (3.1) 

to couple the system to an external gauge field we need an expression for 

G^(x 1 ...x n -x\... x' n , y) = (0 out \T ty\ Xl ) . . . ^(z n tyVi) ■ ■ ■ ^(O^iv)] |0m) , 

(3.2) 

with jv- being the current that couples to the external field. The simplest example of this 
is the quark propagator 

Si(xi - x 2 ) = (0out\T [ipi{x 1 )'ipi(x2)] \0 in ) , (3.3) 
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for a quark of species i, which is gauged to 

S?(x! - x 2 ; y) = (0out\T [M^M^J^y)} \0 in ) • (3.4) 
This also allows to introduce the proper (or amputated) vertex 

Siipf - Pi r = SiipfKipf - pJSifa) , (3.5) 
where we used the definition of Green's functions in momentum spacclike 

G{pi . ..p n ;p'i ■ ■■Pn) = 

J d A Xl . . . d^Xnd^x^ . . . d 4 x' n e- ipiXl . . . e - iPnXn e ip i x 'i . . . e«G(xi . . . x n ; x[ . . . x' n ) , 

(3.6) 

and translation invariance. 

This gauged Green's functions can be obtained by adding a term — j J li {y)A ll {y) and 
taking a functional derivative of the corresponding Green's function 

(3.7) 

A=0 

Therefore, gauging a Green's function is equivalent to taking a derivative and, thus, follows 
the same rules when acting, for instance, on products of functions. 

The coupling of a baryon (through its constituents) to an external field A^ can therefore 
be obtained by gauging the 6-points Green's function G( 6 ). Using equation (1.8) and the 
product rule for gauging a Green's function (and, for the moment, dropping all momentum 
dependence) we obtain 

G^ = G^ - iG^KG - iG K^G - iG KG^ . (3.8) 

We can rewrite this as 

C = (1 + iGoKy 1 (Gfi - iG^KG - iG K^G) 

= G (G^G^Gq 1 - iK") G (3.9) 

where we have used (1.8) to go from the first to the second line. The proper vertex, J M , 
for the coupling of the three-quark system to the field A^ is defined as 

G M = GoJ^Gq . (3.10) 

If the three quarks form a bound state, both before and after interacting with A„, the 
Green's functions G appearing in (3.9) will develop a pole. In this work we are interested 
in virtual* interactions with the external field, in which the initial and final baryons are 



C = 



SA. 



-G 



"This is equivalent to say that the momentum injected by the external field is spacelike. 
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the same. Therefore, taking only the pole contribution, and writing it in terms of the 
scattering matrix T by means of (1.9), Equation (3.9) becomes 

J^ 1 ~ T (Gq — iGoK^Gq) T (3.11) 

or, in terms of the bound-state amplitudes VP (see Equation (2.3)) 

J» = y f (G%-iG K> i G )y i . (3.12) 

We need now explicit expressions for Gq and K^. Using (3.5) and again the product 
rule for gauging the equations, one gets 

Go* = (SiS 2 S 3 r = S»S 2 S 3 + S^Ss + S&Sg = 

s{rlsis 2 s 3 + s 1 s 2 f r%s t 2 s 3 + s 1 s 2 s(r%s l 3 . (3.13) 

To find an expression for K^, one must recall that K is defined as the sum 

K=Y^ K W (3.14) 

£=2,3 

with 

P= Y, K^ S- l . . .S-\ . (3.15) 

cyclic perm. 
£ elements 

Therefore, to gauge K^> one must also gauge the inverse quark propagator S -1 . To do 
this we define 1^ = 0, and thus 

t» = (SS~Y = S^S^ 1 + S^ST 1 )" , (3.16) 

and from here 

(S' 1 )" = -S^S^S- 1 . (3.17) 

We can now write (3.12) as 

= y f ( G% - iGoK i3) ^G -iJ2 ssf sf K i2) '^S i S i + iJ2 sfsfsf K^T^S i S i S i ) 

\ perm. perm. / 

(3.18) 

A diagrammatic representation of this equation is shown in Figure 3.1. 

The expression for J M in the Rainbow-Ladder truncation scheme is considerably sim- 
pler. First of all, the three-body irreducible kernel is absent and so is _fr( 3 )>^. Further- 
more, since the two-body kernel is reduced to a single gluon exchange with all dressings 
embraced in the effective coupling a e //, it does not couple to the external field A^, which 
means that the term i^( 2 )'^ is also absent. This simplified expression is illustrated in Fig- 
ure 3.2. This scheme has been used in [33, 55] to calculate nucleon electromagnetic, axial 
and pseudoscalar form factors using the covariant Faddeev equation. 
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Figure 3.1: Diagrammatic representation of Equation (3.18). 



For spin-3/2 baryons, and restoring all indices and momentum dependence, it reads 



tJq-f 



xi- I I *(y a rz'<y>(pY t ,qY t ,Pf) (S(p{)T^ Pl ,Q)S(p\)) S^(p 2 )S yi (p 3 ) 

p J q 



J 1 } 



a' a 



J p J q 



J p J q 



*aMi\Pi ,Qi , P) - * oL^xKPi >qi >P 



S a > a (pi) {S(pi)T^p 2 , Q)S(j&)) 5y 7 (pa) 

WfrP*. qf \ Pi) ~ * { Xx(P? } A 2 \Pi. 
ScMSpfa) (s(pi)T^ P3 , Q)S(pl] 

<W(pf } , q?\ Pi) ~ ^lx{P?\qf\Pi) 



(3.19) 



where we define 



*S? 7 z = I K Wn'( k ) -W'(P2 - k)Syy,(p 3 + k) * Q( g» 7 »x(p + k,q- k/2, P) , (3.20) 

J k 

as the result of the first term in the Faddeev equation (2.26) and in a similar fashion we 
define ^jd 2 ! and \[d 3 J . "We have introduced the injected momentum Q via the final and 
initial momenta of the interacting quark k 

pf=p K ±% (3.21) 



The relative momenta in the respective terms of (3.19) are, using the definitions in (2.28), 

^ 4 

{2} , Q 

q f/i =^± 4 



{i} ^ /-Q 

{2} /-Q 
pJa } =p±(i-C)| , 



(3.22) 



{3} 
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and since the initial and final states are on-shell, the total momenta are constrained to 
Pf = Pj = —M 2 , with M the mass of the bound state. As is the case for the Faddeev 
equation, we show in Appendix A. 4 that the three terms in (3.19) are formally the same 
when the momentum partitioning parameter is chosen £ = 1/3. 



-9— O — ©- 



J 11 = £ E ^ — 9 — 3 - X - ^ # - 

perm. J „ (_ perm. J ; 



Figure 3.2: Diagrammatic representation of Equation (3.18) in the Rainbow- 
Ladder truncation. 



3.1.1 Normalization of the Bethe-Salpeter amplitudes 

The Bethe-Salpeter equation (2.10), and its Rainbow-Ladder truncated version (2.26), are 
homogeneous and linear integral equations for vl/. Therefore, its solutions are defined up 
to an irrelevant constant. This is not the case, however, for the calculation of the current 
(3.12). Nevertheless, the pole ansatz for the scattering matrix T (see Equation (2.3)) also 
fixes a condition for the normalization of the Bethe-Salpeter amplitudes. This condition 
can be expressed in two different, but equivalent, ways [80, 81]. 
We start with the equation for the scattering matrix 

T = -iK - iKG T . (3.23) 

As explained in Appendix A. 4, the Bethe-Salpeter equation is solved as an eigenvalue 
equation, the real solution corresponding to the eigenvalue A = 1. Thus, we introduce 
a fictitious parameter A, which would correspond the the eigenvalue of the bound-state 
equation when T develops a pole, in the equation above (which we can set to 1 at the end) 

T = -iK - iXKG T , (3.24) 

or 

T = (1 + iXKGoy 1 (-iK) . (3.25) 
Taking a derivative with respect to A 

& = - (i + lX KG )- 2 (iKG Q ) (-iK) = - KG f = TGqT . (3.26) 
ctA (1 + iXKGq) 

Now we want to write this in terms of the Bethe-Salpeter amplitudes. On the bound-state 
pole we write 
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where the dependency on the eigenvalue appears through the bound-state mass M(A). 
Equation (3.26) becomes 

1"J> dM 2 _ ^G ^ 

~(P 2 + M 2 (X)) 2 dX ~ (P 2 + M 2 (A)) 2 ' ( ' 

and from here 

-^ = *G *. (3.29) 

This is the Nakanishi condition [81] for the Bethe-Salpeter amplitudes \&. 

Equation (3.29) can be written in a way such that the eigenvalue A does not appear 
explicitly. Using the Bethe-Salpeter equation 

* = -i\KG V (3.30) 

we obtain, 

i^K' 1 ^ = X^G^ . (3.31) 

Now, since P 2 = —M 2 , the eigenvalue A(M 2 ) can be equivalently considered as a function 
of P 2 . Taking a derivative with respect to P 2 , we get 

d W^&Ai|V ,3,2, 



dP 2 u V dP 2 dP 2 J 

Finally, using the normalization condition (3.29), and afterwards setting A = 1 we find 
another condition for the normalization of the Bethe-Salpeter amplitudes 

' dGp _ . dK- 1 
dP 2 ~ % ^P 2 ' 

This is the Leon-Cutkosky normalization condition [80]. 



•Isa-'-sja-l*-!- ("»> 



3.1.2 Quark-Photon vertex 

To solve Equation (3.19), the quark-photon vertex T^(ki,k 2 ,Q) has to be specified. To 
put the calculation of form factors into the same frame as meson and baryon calculations 
using covariant bound-state equations, we determine the quark-photon vertex as a solution 
of an inhomogeneous Bethe-Salpeter equation [43, 82, 83] 

T ab( k >Q) = Z ^ab + J^ T a»b"( r 'Q) S a"a>(h ~ r)S h n b ,(k 2 + r)K a „ a;b "b(h, k 2 ]r) (3.34) 

with k = (ki + k 2 )/2, Q = k 2 — k\ and Z 2 the quark wave-function renormalization 
constant. We defined 

r"(fci, fc 2 , Q) = -iQ l (2ir) 4 5(4)(k 1 -k 2 + Q)F»(k, Q) (3.35) 

with Qi the electric charge of the quark i. 

The kernel K is the quark-antiquark scattering kernel. Conservation of the electromag- 
netic current (3.19) requires that the bound-state amplitudes are normalized according to 
(3.29) and the the interaction kernel in the quark-photon vertex is truncated in the same 
way [83] as in (3.19). Therefore, to solve the quark-photon vertex we also use a Rainbow- 
Ladder quark-antiquark interaction kernel. 
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Figure 3.3: Diagrammatic representation of the inhomogeneous Bethe- 
Salpeter equation for the quark-photon vertex. 



3.2 Electromagnetic current 



In the last section we derived an expression for the electromagnetic current in terms of 
the photon interaction with the quarks forming a baryon. On the other hand, the form of 
the current is constrained by Lorentz covariance and current conservation to be a linear 
combination of a finite numbers of Lorentz covariants with scalar coefficients (see, e.g., 
[25]). These coefficients are the form factors. 

The electromagnetic current for a spin- 3 /2 particle is characterized by four form factors 
Fi(Q 2 ) [84]. An expression useful for our purposes has been derived in [49]. It reads 



J^(P,Q) = P aa '(P f ) 



+ ((F 3 + W-0^ 



(3.36) 



where P is the Rarita-Schwinger projector (2.6), Pi and Pf are the initial and final baryon 
total momenta, respectively, Q = Pf — Pi is the photon momentum (see Appendix A.l), 
M is the baryon mass and P = (Pf + Pi)/2. The form factors that are measured ex- 
perimentally are the electric monopole (Ge (Q 2 )), magnetic dipole (Gmi(Q 2 ))> electric 
quadrupole (G\e 2 (Q 2 )) and magnetic octupole (Gm 3 {Q 2 )) form factors. They are related 
to the F(s via 

Ge = (l + y ) (^i " rF 2 ) - ~(1 + r) (F 3 - tF a ) , 

G Ml = (l + ^] (F! +F 2 )- 2 ^(1 + t) (F 3 + F A ) , 

V 5 J 5 (3.37) 

G E2 = (F 1 - tF 2 ) - i (1 + r) (F 3 - tF 4 ) , 



G M , = (Fi + F 2 ) - - (1 + r) (F 3 + F A ) , 



with r = Q 2 /AM 2 . It is shown in [84] that if the baryon is spherically symmetric Ge 2 
and Gm 3 must vanish; therefore they measure the deformation of the object. At Q 2 = 
the form factors define the electric charge (e 3 / 2 ), magnetic dipole moment (/U3/2), electric 
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quadrupole moment (Q3/2) and magnetic octupole moment (C3/2) of a spin- 3 /2 particle, 





= G Eo (0), 


^3/2 


= 2^ Gmi(0) ' 


Q3/2 




O3/2 


= 2M 3 G M 3 (0) 



(3.38) 



Once the electromagnetic current is calculated from (3.19), the form factors can be 
extracted using the expressions [49] 

S2 - 2s\ 



Ge 
Gmi 

Ge 2 
Gm 3 



Aiy/T+T ' 

407 (S4 - 2S3) ' 



8ir 2 VTTr 
3i 



2si r + - — rs 2 



16 r 3 



2a 3 r + 



2 

rs4 



(3.39) 



where 



ai = Tr I jwPpupapP^ ? 
S2 = Tr jjwx*j>J j 

s 3 = Tr{j^ 7 £P a P' 3 } , 
s 4 = Tr{j^ Qa 7 ^} , 



(3.40) 



with 7 £ = T£V • 



3.3 Results 

Having solved the covariant Faddeev equation (2.26), we have available not only the baryon 
mass but also all the information about its internal structure through the Faddeev am- 
plitudes In this section we use this information to calculate the electromagnetic form 
factors of the Delta baryon. As in previous chapter, we repeat the calculation for the two 
different models for the effective interaction described in Section 2.4. 

In Figure 3.4 we show our results for the A + form factors. It is worth noting that, due 
to isospin symmetry, the form factors for A ++ , A and A~ differ from those of the A + 
only by a factor corresponding to their charge, as shown in Appendix A. 4. This means, 
in particular, that all form factors for A vanish identically in our approach"!". 

tr This is not the case for spin- 1 /2 baryons. In this case, due to the mixed-symmetry properties of the 
flavor part of the Faddeev amplitude, the neutron has non- vanishing form factors [33] . 
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Figure 3.4: Electromagnetic form factors for the A + for the MT and AFW 
models. The results are compared to lattice data [72] for dynamical Wilson 
fermions at — 384 MeV (green), = 509 MeV (red) and = 691 MeV 
(blue). 
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Due to the short lifetime of the Delta resonance, it is very difficult to study its prop- 
erties experimentally. In fact, the experimental information on the Delta electromagnetic 
properties is restricted to the A ++ and A + magnetic dipoles and, even in these cases, the 
experimental uncertainties make these values unreliable. For this reason, we compare our 
results to a lattice calculation with dynamical Wilson fermions at different pion masses 
[72]. Lattice calculations, however, also suffer from large errors, especially for the electric 
quadrupole form factor. For the magnetic octupole there are no lattice results whatso- 
ever. Moreover, the low photon-momentum regime is inaccessible to lattice calculations. 
An interesting calculation using a constituent spectator quark model has been presented 
in [85, 86]. In this work the authors restrict the Delta to be an admixture of s- and 
d-waves, and fit their parametrization of the wave function to lattice data. Using these 
fits they are able to calculate the form factors, thereby providing a bridge from lattice 
to continuum calculations. When lattice data is not available, we will use these fits to 
compare qualitatively with our results. 

The evolution of the electric monopole form factor Geo with the photon momentum 
Q 2 is shown in the upper-left panel of Figure 3.4. Note that it is ambiguous to compare the 
result of different calculations because each of them gives a different result for the baryon 
mass. Therefore, to eliminate this scale dependence we plot the results with respect to 
Q 2 /M 2 . For Geo, there is a good agreement of our results with lattice data as well 
as a qualitative model independence. As we illustrate in the next section, this result is a 
manifestation of the fact that Geo is relatively insensitive to the detailed internal structure 
of the baryon. 

The calculated charge radius 

<r * 0> - G E0 (0) dQ 2 (3 - 41) 

is shown in Table 3.1 for the MT and AFW models. Compared to the lattice results, 
our values appear considerably overestimated. A possible explanation is the pion-mass 
dependence of the charge radius, which grows as the pion mass approaches the physical 
value from above. Moreover, chiral perturbation theory shows than when the A — > Nit 
decay channel opens the charge radius changes abruptly to a lower value [89] . Since in our 
approach we do not have a mechanism for the Delta to decay, it is therefore reasonable 
that we obtain a higher result for < r EQ > . A combination of these effects could lead to 
a better agreement between lattice results and ours. Incidentally, a number of effective 
models such as Goldstone-boson exchange models and many others (see [86] for a collection 
of results), give also a large value for < r 2 EQ >. 

In the upper-right panel of Figure 3.4 we plot the magnetic dipole form factor. The 
dimensionless magnetic dipole moment Gmi(O) is 2.22 and 2.33 for MT and AFW, 
respectively. This is to be compared with the experimental value 3.54^^2 [14, 88] or 
the lattice extrapolated result 2.35 ± 0.16 [72] (see Table 3.1). For the Q 2 -evolution 
we observe again a qualitative model independence, but the behavior in both cases is 
significantly different to that of lattice data. 

The deformation of the Delta is signaled by the electric quadrupole and magnetic oc- 
tupole form factors. In the non-relativistic limit, a negative or positive electric quadrupole 
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F-MT F-AFW 


DW1 DW2 DW3 


Exp. 


< r 2 E0 > (fm 2 ) 


0.67 0.60 


0.373 (21) 0.353 (12) 0.279 (6) 




Gmi(O) 


2.22 2.33 


2.35 (16) 2.68 (13) 2.589 (78) 


3-541^ 



Table 3.1: Comparison of results for the charge radius < r% > and for 
Gmi(0) (oc fi). We compare our results for the MT model (F-MT) and for the 
AFW model (F-AFW) to a lattice calculation with dynamical Wilson fermions 
at = 384 MeV (DW1), m w = 509 MeV (DW2) and m n = 691 MeV (DW3) 
[72, 87]. For Gmi(0) we also compare to the experimental value [14, 88]. 



moment corresponds to an oblate or a prolate charge distribution, respectively. In our cal- 
culation, and for both effective interactions, Ge2 starts with a positive value at low-Q 2 , 
changes sign and reaches a minimum at Q 2 jM\ ~ 0.1 and then approaches zero from the 
negative sign, in agreement with lattice data. The behavior at low-Q 2 is in contradiction 
with the constituent-quark calculation in [85, 86], with an extrapolation of lattice data 
[72] and with a covariant quark-diquark calculation [49]. In this respect it is important 
to point out the 1/Q 4 factor in the expression for Ge2 in (3.39). This entails that very 
precise cancellations must take place in the denominator in (3.39) to cancel the Q 4 factor 
and get finite results close to Q = 0. This is difficult to achieve numerically and, in fact, we 
needed to solve the Faddeev equation at the highest precision available to obtain reliable 
results at a reasonable small Q. Using as a measure of the numerical accuracy the value 
of the imaginary parts of the form factors (which, ideally, should be zero), we plot those 
values for which the ratio between the imaginary and the real parts of all form factors is 
smaller than 10 -9 ; the minimum Q 2 -value we are able to achieve in this way is ~ 0.01. In 
any case, this low-Q 2 behavior must be considered with caution and further studies are 
required. If this feature of Ge2 is confirmed, it will be a true prediction relying on the 
high orbital angular momentum components of the Faddeev amplitude, as shown in next 
section. 

A similar situation appears for the magnetic octupole form factor Gm3- In this case 
we obtain a small, but non-zero, positive form factor at high-Q 2 , with a similar behavior 
for both models. At low-Q 2 the models differ, but both seem to feature a pump and a 
zero crossing (i.e. a change of sign of the form factor). At Q ~ 0, a factor l/Q 6 must be 
canceled by the numerator in (3.39) and, therefore, the result here must be taken with due 
caution. It is remarkable, nevertheless, that this behavior of Gm3 is qualitatively similar 
to that of [85, 86]. These features at low-Q 2 are also absent if one ignores the subleading 
components in the Faddeev amplitude (see next section). 

In summary, we find a reasonable model independence of our results, which gives us 
confidence that we can make statements about the Rainbow-Ladder truncation of the 
baryon BSE. The models slightly disagree at low-Q 2 for the electric quadrupole and the 
magnetic octupole form factors but, as we discussed above, it is not clear at the moment 
to what extent this is a numerical artifact. 

To conclude, it is noteworthy that in a covariant approach there is no freedom to choose 
the relative importance of the different quark-spin ar quark orbital angular momentum 



40 



CHAPTER 3. SPIN-3/2 ELECTROMAGNETIC FORM FACTORS 



components. Once the parameters of the interaction are chosen, the contribution of each 
component is determined by the Faddeev equation. This is in contrast with, for exam- 
ple, the constituent-quark calculation in [85, 86] where the different angular-momentum 
contributions have to be fixed a priori and fitted afterwards to lattice data. 

3.3.1 Role of the subleading components 

The Poincare-covariant description of spin- 3 /2 baryons forces us to use a basis of dimension 
128 for the Faddeev amplitudes-'", as shown in Appendix A. 2. This basis can, in turn, be 
classified in terms of the quark spin (s) and orbital angular momentum {£). At first sight, 
this seems an exceedingly complicated way of describing a baryon. In fact, for spin-3/2 
baryons the s-wave sector (s = and £ = 0) is dominant. To illustrate this we have plotted 
in Figure 3.5 the dominant amplitudes in each of the (s,£) sectors (specifically, we plot 
the dominant Chebyshev moment, as defined in (A. 57), of this amplitude) for the Delta 
in the MT model, although a the same reasoning applies for the AFW model. 

To evaluate the relevance of the subleading components we performed the following 
calculation. After solving the Faddeev equation for the MT model and using all 128 basis 

(i) S 

elements we isolate those components fi^ m such that s 



with {pi,qj} the quadrature points used in the numerical integration and N po i nts the 
number of quadrature points used. Note that this condition especially implies that each 
of these components is a factor of, at least, 1CP 7 smaller than the dominant one. 

We then remove this components from the basis and solve again the Faddeev equation 
using the remaining basis subset. The resulting Delta mass is 1.30GeV, to be compared 
with the result of the full calculation Ma = \.22GeV. That is, the difference between 
both results is about 6%, and the mass with the reduced basis is only 2.3% above the 
physical Delta mass. This might suggest that, indeed, only the dominant components (i.e. 
s-waves) are necessary in the calculation (especially if, as we concluded in the last chapter, 
the results in Rainbow-Ladder are expected to be precise up to ~ 10%). 

The result of the calculation of the Delta electromagnetic form factors using the reduced 
basis is shown in Figure 3.6. The result for the electric monopole Geo and magnetic dipole 
Gmi form factors are similar to the ones in the full calculation. The electric quadrupole 
form factor Ge2 is also qualitatively similar for the full and the reduced calculation, 
although a bit featureless in the latter case. Finally, the behavior of the magnetic octupole 
form factor Gm3, which is indicative of the baryon deformation, is very different to the 
corresponding behavior in the full solution. Therefore, we conclude that although a rough 
description of baryon properties is achievable using a simplified setup, the study of fine 
details such as deformation requires to keep all components, as dictated by symmetry. 

"'"For spin- 1 /2 baryons the dimension of the basis is 64. 

§ Note that the normalization of the Faddeev amplitudes for this exercise is arbitrarily fixed such that 
foooiP 2 ~ 0,q 2 = 0) = 1- This normalization is not related to the physical normalization given by (3.29), 
which is used for the calculation of form factors. 




(3.42) 
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Figure 3.5: Dominant Chebyshev moments in each of the (s, L) sectors (we 
change here the notation for the orbital angular momentum from I to L for 
clarity), as a function of p and for different values of q in the Maris- Tandy 
model. The amplitudes are normalized such that /ooo(<Z 2 = 0,p 2 =0) = 1. 
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Figure 3.6: Result for the Delta form factors when the subleading components 
are removed. 
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Yet another manifestation of the importance of subleading components comes from 
the comparison of the results presented here with those obtained using the quark-diquark 
approximation [49] . The description of baryons using diquark effective degrees of freedom 
naturally misses some of the elements of a full three-body description. As a consequence, 
the qualitative behavior of Ge2 and Gm3 in the quark-diquark and in the Faddeev ap- 
proach is different. 
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Chapter 4 

A note on glueballs 



The quark model allows two families of hadrons: mesons, formed by a quark and an 
antiquark and baryons, formed by three quarks. States not describable by the quark 
model are called exotic states. In particular, a longstanding prediction of QCD is that 
there should exist bound states formed by gluons only, so-called glueballs [36]. 

An enormous experimental effort has been put on studying meson spectroscopy (see, 
e.g., [90] for a review) resulting on the discovery of exotic hadrons at Belle [91, 92] and 
Babar [93] collaborations, but for the moment there is no direct evidence of glueballs. 
The reason for this is that the nature of these exotic states is not well understood; one 
possibility is that they are normal quark-antiquark states mixed with glueballs since they 
have the same quantum numbers. Glueballs could also be hidden in the background of 
other mesons which are more abundantly produced. To improve this situation, there are 
several planned experimental facilities with the glueball search as one of the primary goals. 
GlueX at Jefferson Laboratory (USA) is expected to run in 2014 and will focus on mapping 
the spectrum of exotic states and on the search of light-quark hybrids which may serve 
to study glueballs via their decays. PANDA at FAIR/GSI (Germany) will study non- 
exotic hadrons up to charmonium states and from their decays will search gluonic states. 
However, a deeper theoretical understanding of glueball and other exotic-state formation 
must proceed in parallel with these experimental endeavors. For example, a calculation of 
glueball masses from first principles and the description of possible mixtures of glueballs 
and standard quark-antiquark bound states would be desirable. 

Glueballs are also interesting from a purely theoretical point of view. The reason why 
QCD predicts the existence of glueballs is because it is a non-abelian gauge (or Yang- 
Mills) theory, which in particular implies that, as opposed to abelian gauge theories like 
Quantum Electrodynamics, the gluons (or gauge bosons) can interact among themselves 
and therefore can form pure gluonic (i.e., without quarks) bound states. Thus, glueballs are 
an excellent probe into the non-perturbative, or low-energy, regime of non-abelian gauge 
theories. However, most of the glueball spectrum calculations have been done by modeling 
QCD with effective degrees of freedom (see, e.g., [94] for a review). The first calculations 
of glueball masses were performed using bag models with constituent massless gluons 
[95-100], where gluons are confined into a potential bag. Models with massive constituent 



45 



46 



CHAPTER 4. A NOTE ON GLUEBALLS 



gluons have also been studied [101-104], however they introduce spurious states due to the 
unphysical longitudinal degree of freedom of a spin-1 constituent gluon. Other approaches 
used to calculate the glueball spectrum are QCD sum rules [105-107], flux-tubes [108] or 
AdS/QCD [109]. There is no general agreement about the glueball masses among different 
approaches, although some qualitative features can be extracted [110, 111]. Moreover, it 
is not clear whether these models capture all features, or any, of QCD dynamics. 

A connection of the constituent gluon models with QCD was attempted in [112] by 
deriving an effective Hamiltonian from the QCD Lagrangian. Later, a calculation of the 
low-lying glueball masses based on continuum QCD was presented in [113]. Specifically, 
the authors construct a Fock space of constituent quarks and gluons and an instantaneous 
potential between them using Coulomb gauge QCD. So far, only lattice QCD calculations 
provide values of the glueball spectrum using the fundamental degrees of freedom of QCD 
[114-118]. These calculations have nevertheless some drawbacks. First of all, most of them 
are performed in the so-called quenched approximation, in which the quantum fluctuations 
involving quarks are neglected, and there is no consensus about the effect of including those 
quark loops [114, 117, 118]. Also, common to all lattice calculations, are the problems 
of discretization and finite-volume effects which indeed yield to different glueball masses 
in different calculations. Finally, it is difficult from lattice calculations to unravel how 
gluonic bound-states are formed. 

To improve the understanding of how glueballs emerge from QCD it is convenient to 
have an approach complementary to lattice-QCD. This could help to identify the leading 
mechanisms responsible for the formation of bound-states out of the elementary degrees 
of freedom of the theory and therefore simplify the theoretical treatment by keeping only 
those mechanisms in the calculations. This is especially important if one eventually at- 
tempts to study the mixing of gluonic and ordinary bound-states. This complementary 
approach is provided by functional methods like Dyson-Schwinger or Functional Renor- 
malization Group equations. 

However, so far the calculation of glueball masses using a covariant Bethe-Salpeter 
approach has been very little explored. In this chapter we propose a BSE for the simplest 
glueball, namely with no admixture of quarkionic states and formed only by two valence 
gluons. 

4.1 Glueball Bethe-Salpeter equation 

Throughout this work we have stressed that the quark-antiquark or the three-quark 
Green's functions develop a pole when the system forms a bound state. We have used this 
fact to go from the Dyson equation for the Green's function (1.8), or for the scattering 
matrix (1.10), to the Bethe-Salpeter equation for the bound-state amplitude. These equa- 
tions rely upon an, in general, unknown interaction kernel K which, nevertheless, can be 
expanded diagrammatically in a (more or less) straightforward way. 

Green's functions are equivalently defined by Dyson-Schwinger equations. Thereby, 
the equation is completely fixed, without the presence of undetermined interaction kernels. 
Instead, the equation depends on the knowledge of higher order Green's functions and, 
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therefore, one has to solve an infinite and coupled system of Dyson-Schwinger equations. 
In any case, a truncation (or ansatz) must be chosen either for the interaction kernel K 
or for the high-order Green's functions. 

In this section we use the pole assumption to derive a Bethe-Salpeter equation for 
a system of two gluons forming a glueball, using the Dyson-Schwinger equation for the 
two-gluon four-points Green's function. Before that, however, we illustrate the procedure 
using the meson Dyson-Schwinger and Bethe-Salpeter equations as an example. 

4.1.1 Warm-up: meson Bethe-Salpeter equation 

The evolution of a quark-antiquark system is described by the four-points Green's function 
G q q(x' 1 ,x' 2 ,x 1 ,x 2 ) = (0\q(x' 1 )q(x' 2 )q(x 1 )q{x 2 )\0) , (4.1) 
or, more conveniently for our purposes, by its momentum-space counterpart 

Gqq(p' 1 ,P 2 ,Pl,P2) ■ (4.2) 

The derivation of the Dyson-Schwinger equation for this Green's functions is very involved 
and for these reason we use DoFun, a Mathematica™ package described in [119, 120]. 
The result (more precisely, the equation for the amputated Green's function, or scattering 
matrix T) is shown in Figure 4.1. 




Figure 4.1: Dyson-Schwinger equation for the quark-antiquark four-points 
amputated Green's function. Big blobs denote full Green's functions. 

To obtain a Bethe-Salpeter equation we proceed as follows. Whenever the incom- 
ing quark and antiquark momenta pi and p2 are such that (pi + P2) 2 = P 2 = — M 2 , 
where M is some meson (i.e. bound state) mass, the momentum-space Green's function 
Gqq(p'i,p' 2 ,pi,P2) will have a pole. Of the diagrams in Figure 4.1, only those in which 
the incoming quark and antiquark enter a quark-antiquark four-points Green's function* 

*In principle, a skeleton expansion (that is, in terms of primitive Green's functions) of higher-order 
Green's function may contain sub-structures which consist of quark-antiquark four-points Green's functions 
and should, therefore, be included as well. We ignore here this possibility. 
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in the s-channel will potentially develop a pole when the bound-state is formed. This 
procedure has been applied in [121-123] to study an explicit implementation of the BRST 
quartet mechanism [124]. 

It is clear that only the first diagram in Figure 4.1 fulfills the conditions stated above. 
Keeping only this diagram and expanding, as usual, the amputated Green's function 
around the bound-state pole to introduce the Bethe-Salpeter amplitudes \I> 



T 



P 2 + M 2 ' 

we obtain the Bethe-Salpeter equation for the meson, as shown in Figure 4.2. 



(4.3) 



Figure 4.2: Bethe-Salpeter equation for the meson, 
propagators or vertices. 



Big blobs denote full 



The conventional Bethe-Salpeter equation for a meson would be 

V = K q jGf ) * (4.4) 

with an undetermined kernel K. The relation between this equation and the one derived 
here (Figure 4.2) should be obtained, in principle, by using a skeleton expansion of the 
full quark-gluon vertex. It is not clear, however, how this expansion would reproduce all 
diagrams in K'. It is obvious, nevertheless, that within the Rainbow-Ladder truncation 
scheme both equations are identical. 



4.1.2 Glueball Bethe-Salpeter equation 

As illustrated in previous section, the first step to study a given state in a quantum field 
theory is to determine the operator which, acting on the Hilbert-space vacuum, creates 
that state. In the case of QCD, observable particles correspond to gauge invariant local 
operators. Based on this idea, an analysis of the possible glueball quantum numbers was 
performed in Refs. [36, 111]. It is also possible to mimic the simplistic picture of the quark 
model and classify the lightest glueballs as composed of two or three gluons. This widely 
used classification makes implicitly a difference between the so-called valence gluons and 
the gluons which bind those into a bound-state. In this project we assume this picture 
and will concentrate on two-gluon glueballs. 

The simplest operator describing a glueball is: 

(0\Fr(y)F^(x)\0) . (4.5) 

^Possibly including some of the terms in a skeleton expansion of higher-order Green's functions in the 
equation. 
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This operator provides information about the gluon propagation as well as its self- interactions. 
In particular it contains a four-gluon operator which will give rise to the glueball. 

We use DoFun to calculate the Dyson-Schwinger equation for the two-gluon four-points 
Green's function. The output contains 79 diagrams which we will not show here. Applying 
the procedure discussed in previous section we end up with the Bethe-Salpeter equation 
shown in Figure 4.3, which consists of five diagrams. 




Figure 4.3: Proposed Bethe-Salpeter equation for the glueball. Big blobs 
denote full propagators and vertices. 

The glueball Bethe-Salpeter amplitude can be written, in general, as 

< lW ... M >,^) , (4-6) 

where the [i, v Lorentz indices and the a, b color indices refer to the valence gluons and 
the rest of Lorentz indices describe the angular momentum of the glueball. The relative 
and total momenta are 

P = Pi + P2 = Pi + p 2 , 

Pi —P2 (4-7) 



with pi and p 2 (and p\ and p' 2 ) are the valence-gluon momenta. 

A remarkable, and somewhat unexpected, feature of this equation is that ghosts do not 
appear explicitly. There are essentially two ways ghost lines could be explicit. The first 
one is through a skeleton expansion of higher order Green's function in the four-gluon DSE 
but, as explained before, our truncation scheme enforces to neglect these diagrams. The 
second possibility is to assume a priori that ghosts must appear explicitly and study the 
four-gluon and two-gluon - two-ghost coupled system. It is important to notice, however, 
that in the equation we propose ghosts can play an important role in the dynamics of the 
system. For example ghost propagators and ghost-gluon vertices are essential elements in 
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the gluon DSE (see Figure 4.4). The same can be said about the DSE for the three-gluon 
vertex. 




Figure 4.4: Gluon DSE. Wiggly, solid and dotted lines represent gluon, quark 
and ghost propagators, respectively. Dressed Green's functions (propagators 
or vertices) are denoted by blobs. To solve the equation one needs to know the 
dressed quark and ghost propagators and three-gluon, four-gluon, ghost-gluon 
and quark-gluon vertices which are given by the corresponding DSEs that in 
turn depend on other Green's functions. 

The Bethe-Salpeter equation obtained in this way is extremely complicated. First 
of all, the diagrams with two gluon-loops will be UV-divergent and some regularization 
procedure must be used. On the other hand, the third diagram in Figure 4.3 contains the 
off-shell full four-gluon vertex. Since the BSE is derived from the full four-gluon vertex 
DSE assuming this vertex shows a pole, to include this diagram we need an ansatz which 
is consistent with the pole assumption. This seems a highly intractable problem. 

In the absence of a way to deal with those problems, we simply ignore the abovemen- 
tioned diagrams and define our glueball BSE with only the first two diagrams in Figure 
4.3. However, we still face another problem. To see this we expand the amplitude (4.6) in 
a basis 

s:ab 

with £\ /(0 T W the spin part of the amplitude and the color part , with {a, b} 

v N c 1 

indices in the adjoint representation, is fixed by requiring the glueball to be a color singlet. 
Since gluons are bosons, the total amplitude (4.6) must be symmetric upon the interchange 
of the valence-gluon indices and momenta. Therefore, the color part 5 ab forces the spin 
part to be symmetric. Now, consider one has a symmetric amplitude as the input for 
the second diagram in Figure 4.3. Since only one of the three-gluon vertices is full, the 
output of this diagram will not be symmetric and therefore this term cannot appear in the 
glueball BSE if we want to keep Bose-symmetry. We propose to fix this problem by using 



4.2. STEPS TO SOLVE THE GLUEBALL BETHE-SALPETER EQUATION 



51 



full vertices for both three-gluon vertices in this diagram. The resulting glueball BSE is 
shown in Figure 4.5. 




Figure 4.5: Simplified Bethe-Salpeter equation for the glueball. Big blobs 
denote full propagators and vertices. 

A realistic description of glueballs must deal with the mixing of them with quark- 
antiquark states. This is an issue of exceptional importance in the experimental search 
of glueballs and so far there is no clear theoretical account for it. In the DSE/BSE 
formalism the mixing of quarkionic and gluonic states would be studied by solving a 
coupled system of BSEs. The equation derived here assumed that only the four-gluon 
Green's functions develops a pole. However, in the case of a mixed quarkionic-gluonic 
state, the poles corresponding to a bound-state would appear simultaneously in the four- 
gluon, four-quark, and the quark-quark-gluon-gluon Green's functions. Therefore the BSE 
should be modified. 

It is also important to remark here that the full Green's functions contain information 
about bound-states with all possible quantum numbers. A specific choice of quantum 
numbers for the bound-state of interest is done via the symmetries and structure of the 
BS amplitudes. 



4.2 Steps to solve the glueball Bethe-Salpeter equation 

To write an explicit expression for the glueball Bethe-Salpeter equation, we need to in- 
troduce the definitions of the different Green's functions involved. The tree-level gluon 
propagator in Landau gauge is 



with Zj, the gluon wave- function renormalization constant and Z(p 2 ) the gluon propagator 
dressing function. The renormalized four-gluon vertex is 

+ f ade f bce (g^g pa -g w gua) 

(4.10) 

where f abc are the SU(N C ) structure constants, Z4 the four-gluon renormalization constant 
and all momenta are considered as outgoing. Finally, for the renormalized full three-gluon, 
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we assume that the color structure is the same as the tree-level one and we write it as 

^ C p (k,P,q) = -i 9 Z 1 (27r) 4 S^(p + q + k)T pil/p (k,p,q) (4.11) 

so that for the bare vertex one has T^ up (k,p,q) = T$ p (k,p,q) = {k-p) p g pu + (p-q) p g up + 
(q — k) v g pp . Z\ is the three-gluon renormalization constant and, again, all momenta are 
taken as outgoing. 

It is also necessary to define a covariant basis t^' for the expansion (4.8). A systematic 
derivation of such a basis for glueballs of any spin has been described in [125]. We shortly 
describe here the main ideas to obtain this basis. The first step is to couple the two 
valence gluons (that is, their Lorentz indices) into a Lorentz invariant object. The only 
possibilities are 

g »v ^ £Ppv € PV»v . ( 4-12 ) 

The scalar glueball ++ is represented by g^ v /2 whereas the appropriate basis element for 
the pseudoscalar glueball CT + is 



M-\f> ■ P) pf>P°e pr7 ^ , M = yj{p ■ P f (l - (p ■ P) 2 ) , 



(4.13) 



where the factor (p ■ P) is included to obtain a state of positive charge conjugation, since 
this operation amounts to a sign flip in the relative momentum for a system of two gluons. 
Basis elements for higher angular-momentum states are constructed with tensor products 
of the aforementioned Lorentz scalars with Lorentz tensors with the appropriate number 
of degrees of freedom. To represent a state of angular momentum J, the tensor must have 
2 J + 1 independent components. As described in [125], this is achieved with tensors that 
are symmetric, traceless in all pair of indices and orthogonal to the total momentum in 
all indices. The construction of these states grows very fast in complexity. 
With these definitions, the first diagram reads 

1 ^2 f d 4 p' _ m , 6 ab 



x [f abe f cde (g^-g^x) 
+ f ace f bde (g^gxa-g^gux) 

+ f ade f bce {g^g^ - g^g^ 



x^, m ...,. 7 (p'^)^ f /°' ) (p' 2 ,P'-P;P 2 ) 



(4.14) 



where we used the Slavnov-Taylor identity Z4 = Z%Zg and introduced the renormahzed 
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coupling constant g(fi) = Z g g. The color traces for the three terms in the integral are 

^ab jabe jcde gee' ^dd' jc' a" jab jabe jcce q 



S ab^aeejbde d ee' § dd' 6 e'd' = § ab jaee jbee = gobgab^ = ^2 _ ^ ^ ^ ^ 

gab jade jbee gee' gdd' gc' ' d' gab jace jbee gabgabj^ (./V^ l)iV 

where we used the identity f acd f bcd = N c 5 ab . 
The gluon-exchange diagram reads 



x f brs r^( Pl , -p[, k)r tm T^ p {v2, —k, -p' 2 )Df s (k)D^p[)D^(p' 2 ) 

gul 



x 4i], 1 ...^p',P)^=f^( P ' 2 ,p' ■ P;P 2 ) , 



(4.16) 



where we used the Slavnov- Taylor identity Z\ = Z g Z^ 2 and g(fi) = Z g g. The color trace 
gives 

gab jatm jbrsgstgmlgruglu gab jatm jbmt gab jatra jbrat gab gab (_/V^ l)iV 

('1-17) 

It is important to notice that, if one is given with expressions for the propagator and 
vertex dressing functions, Equations (4.14) and (4.16) have only g(fi) as a free parameter. 

For a self-consistent glueball calculation in the DSE/BSE formalism it is essential to 
solve the gluon propagator from the gluon DSE (see Fig. 4.4). The gluon propagator has 
been studied numerically in [62, 63, 126] using several approximations. In both works it 
was concluded that the gluon propagator in the complex plane has a non-trivial analytic 
structure, although in contrast to the quark propagator used in previous chapters, it has 
only branch-cuts with absence of poles. Non-analiticities of the gluon propagator will 
presumably be of relevance in the resolution of the glueball BSE. 

To solve the gluon DSE one needs to know the dressed quark and ghost propagators and 
three-gluon, four-gluon, ghost-gluon and quark-gluon vertices. Exact expressions for these 
quarks and vertices are generally unknown and models, at least for the vertices, must be 
used. Restrictions on those models can be based on infrared analysis of Yang-Mills vertices 
[127, 128] and on lattice studies [129, 130], as well as from the underlying symmetries 
of QCD (formalized by the Slavnov- Taylor identities). A self-consistent solution of the 
glueball BSE is, nevertheless, beyond the scope of this work. 
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Summary 



The calculation of hadron properties from QCD remains an open issue due to the fun- 
damental non-perturbative nature of the problem. Covariant Bethe-Salpeter equations, 
in combination with Dyson-Schwinger equations for QCD Green's functions, constitute 
an excellent tool for hadron studies in continuum quantum field theory. They allow, in 
particular, a systematic analysis of the mechanisms in QCD that lead to the formation of 
bound states. Moreover its applicability is not restricted, in principle, to any momentum 
or quark- mass range, which is a limitation in other approaches such as lattice QCD. 

The limitations of the approach come from the fact that a complete solution of the 
Bethe-Salpeter equation would require to solve an infinite set of coupled Dyson-Schwinger 
equations for the QCD Green's functions. A realistic calculation thus requires a truncation 
of the system, which means that of all the possible correlations among particles provided 
by QCD, one considers only a subset. Any truncation, in turn, induces the necessity 
of some modeling. It is therefore necessary to disentangle from the results the model- 
dependent and the truncation-dependent features. The latter carry the information about 
the relevance of the interaction terms considered in the calculation. 

In this thesis we focused on the study of baryon properties using a three-body covariant 
Bethe-Salpeter calculation. We used the Rainbow-Ladder truncation for the quark-quark 
kernel. This truncation consists of a vector-vector dressed-gluon exchange between quarks. 
Since we do not solve the corresponding Dyson-Schwinger equations for the gluon prop- 
agator and the quark-gluon vertex, a choice of a model for their dressing functions is 
necessary. Moreover, we neglect three-body irreducible interactions. 

We performed our calculations using two different models. One is the Maris- Tandy 
model, which is a purely phenomenological one, designed to provide dynamical chiral 
symmetry breaking and a good description of ground-state meson properties. The second 
model, which we called the Alkofer-Fischer- Williams model, attempts to capture some of 
the QCD dynamics in the infrared and, therefore, might be a closer approximation to a 
calculation from first principles in QCD. 

Using those two models, in Chapter 2 we calculate the masses of spin-1/2 and spin- 
3/2 baryons for a range of quark masses from u/d up to b quarks. The conclusion of 
this chapter is no strong model dependence is manifested in the baryon spectrum. The 



55 



56 



CHAPTER 5. SUMMARY 



difference between models, and between the models and the physical or lattice values is 
always smaller than 10%. Our results also indicate that irreducible three-body interactions 
are subdominant in baryons. However, a more precise determination of its relevance 
would require to study also quark-quark interaction terms beyond the Rainbow-Ladder 
truncation. 

In Chapter 3 we study the electromagnetic form factors of the Delta(1232), again 
using the two models for the effective interaction. In this context it is important to 
remark that, in a covariant approach, all partial waves allowed by symmetry are present 
in the calculation and that their relative importance is dictated by the dynamics. In the 
case of spin- 3 /2 baryons, this implies the presence of s-, p-, d- and f-waves. The most 
important result in this chapter is that the Delta has non-zero electric quadrupole and 
magnetic octupole form factors, which is a signal of deformation from sphericity. Although 
the Delta(1232) is dominated by s- waves, we also show the importance of the subleading 
components in the electromagnetic properties. In the calculated form factors we also 
observe a qualitative model independence. 

Most of the experimental information about the electromagnetic properties of the Delta 
and the details of its deformed shape is obtained, however, from the electromagnetic decay 
N —7- A7. In a covariant approach, the nucleon and Delta Bethe-Salpeter amplitudes and 
the quark-photon vertex used for the form factor calculations are also the elements needed 
to study this process, which will be the object of future work. 

Covariant Bethe-Salpeter equations are not limited to the calculation of hadron prop- 
erties in QCD. They can be extended to study bound-states in any quantum gauge-field 
theory and also to the study of exotic states in QCD. In Chapter 4 we take a first step in 
the calculation of glueball properties in the BSE/DSE framework. We propose a Bethe- 
Salpeter equation for glueballs based on the four-gluon Dyson-Schwinger equation in QCD. 
However, an explicit solution of this equation is very difficult and is part of current and 
future investigation. 
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A.l Conventions and reference frames 



In this thesis we work in Euclidean spacetime. Four vectors can be expressed using hy- 
perspherical coordinates: 



P^ = VP 2 



/ sin a sin (3 sin ip \ 
sin a sin (3 cos <p 
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I \J\ — z 2 \J\ — y 2 simp \ 
\J\ — z 2 a/1 — y 2 cosy? 
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v 7 ! - z 2 y 



(A.l) 
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in which the four-momentum integration is written as 

In the resolution of the Faddeev equation one has two external (i.e. they are not 
integration variables) momenta p, q. In this a case one can choose a reference frame such 
that they can be written as 
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(A.3) 



With this choice and on the baryon's rest frame = (0, 0, 0, i M), the normalized mo- 
menta pr^ and q t ^ needed for the construction of the basis for the Faddeev amplitudes 
(see AppendixA.2) simplify to 
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(A.4) 



For the calculation of form factors, however, the baryon's rest frame is not convenient. 
We use instead the Breit frame (or z-Breit frame) . In this reference frame the baryon has 
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initially a four-momentum Pj* and interacts elastically with a photon with four-momentum 
Q to give a baryon with final four- momentum Pj, such that 



if = -Pjf , = P» - Jf , 








\ 



±i Q/2 M 



/ \ 



IQI 
V o / 



(A.5) 
(A.6) 



Note that Pf = P 2 = — M 2 as it should be, since the baryon is on-shell in both cases. 



/ 

For the Euclidean Dirac matrices we use the convention 



7 



1 
1 



r 





i a 1 



-i a 




7 



1 
-1 



7 



1 

1 



with o l the Pauli matrices. This choice of matrices is hermitian 7 M = (7^) 



(A.7) 



A. 2 Covariant decomposition of spin- 3 /2 baryon amplitudes 

The Faddeev amplitudes are tensors that describe the baryon in terms of the valence 
quarks 

* ~ (01^^410) • (A.8) 
where is a baryon creation operator. Schematically they can be decomposed as 

SPIN ® FLAVOR ® COLOR . 

The color part of the amplitude is fixed by the fact that the quarks must combine into a 
color singlet, and therefore for a system of three quarks it is given by the antisymmetric 
tensor 

COLOR = e -^B- (A.9) 

being A, B and C the quark color indices (in the fundamental representation) and the a/6 
is included to normalize it. As for the flavor part, it is given by the usual quark- model 
flavor states. 

The structure of the spin part depends on the baryon of interest. In general it can be 
written as 

^a^x(p,q;P) (A. 10) 

where p, q and P are the relative and total momenta, as defined in (2.28), 

P 1 - C 

p = (1 - C)P3 - C(pi +P2) , pi =-q--^-\ — y~ P ' 

P2-P1 P 1 - C T3 (A. 11) 

Q , P2 =l- 2 + ^ r P> 

P = pi + p 2 + p 3 , p 3 = v + C i P, 
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with pi, P2 and p% the quark momenta and £ the momentum partitioning parameter. 
The indices a, f3 and 7 represent the valence quark and the generic index I indicates the 
baryon of interest. For spin- 1 /^ baryons this index refers to a Dirac field and, therefore, 
the Faddeev amplitude (since the color and flavor parts are trivial, when we talk about 
Faddeev amplitudes we refer to its spin part, unless otherwise stated) is a rank-4 Dirac 
tensor. For spin-3/2 particles the index I refers to a Rarita-Schwinger field and therefore 
the Faddeev amplitude is a mixed tensor with four Dirac indices and one Lorentz index. 

These tensors can be conveniently expressed in terms of a basis (see Equation (2.31)). 
It can be shown that the positive parity and positive energy subspace of the spin-!/2 basis 
contains 64 linearly independent elements [50]. The positive parity and positive energy 
subspace of the spin-3/2 basis contains 128 elements [34]. In what follows, we show the 
steps for the construction of a covariant basis for the Faddeev amplitudes. 

The construction of the corresponding basis is in principle straightforward, but very 
cumbersome due to the large number of indices and basis elements involved. For this 
reason is most convenient to perform all the calculations using some symbolic programming 
language like Wolfram's Mathematica™ . The most significant aspects of this construction 
are: 

• It is independent of any approximation in the three-body Bethe-Salpeter equation. 

• Only Poincare covariance as well as parity invariance are needed to construct the 



• The basis includes all possible internal quark-spin and orbital angular momentum 
values that can lead to the final spin of the baryon of interest. 

The last point is very important. In a complete covariant calculation there is no freedom to 
choose the spin and orbital angular momentum composition of the baryon. All possibilities 
must be included, this is dictated by symmetry, and the relative importance of each of 
them will be determined by the interaction. 

A basis for a rank-four Dirac tensor can be constructed out of the following linearly 
dependent elements, 



basis. 




/ 



(r^r.o^P) 



(A.12) 



where z, j = 1 ... 4 and 




(A.13) 



where the hat denotes a normalized vector 



(A.14) 
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and, for convenience, the unit transverse four-vectors pr and q\ are used to contract the 
Poincare indices 



with 



the transverse projector. We also defined 

n a (P) = A a (P) 7 5 C ® A+(P) 



(A.15) 
(A.16) 

(A.17) 
(A.18) 



with a = ±, C = 7 4 7 2 is the charge conjugation- matrix, A ± (P) = (1 ± ^) /2 is the 
positive- and negative-energy projector. These are 128 elements of which one can check 
that only 64 are linearly independent. Its structure and partial- wave decomposition can 
be found in [50]. 

Using the elements in (A. 12) we can construct the building blocks of a mixed Dirac- 
Poincare basis in the following way 



M 9 - 



Mf 3 



M q - 
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( It® 1 \ 



(M§)(10P M1/ ), 



(A.19) 



where M G {S,P,V, A} are the basis elements defined in (A. 12) and P 1 *" is the Rarita- 
Schwinger projector for positive-energy particles 

Pf (P) = A+(P) (t^ - l -^ TlT ^ , (A.20) 

with 7^ = Tp U j u . The set (A.19) contains 384 elements, but it can be checked that 
only 128 of them are linearly independent. Of course it is physically irrelevant which 128 
elements are chosen to form a basis. Nevertheless, from a naive quark-model point of view, 
the s-wave components (i.e. the relative-momentum independent basis elements, such as 
Sf : ) will play a dominant role and therefore is convenient to include them explicitly in the 
basis. 

To this end, the basis elements can be classified with respect to their quark-spin and 
relative orbital angular momentum in the baryon's rest frame, as explained in the following 
subsection and summarized in tables A.l and A. 2. For the purpose of this classification, 
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we found the following choice of linearly independent elements convenient: 
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(A.21) 



qP q9 
S 3i > S 4j 
pP p<? 
*3j ' r 4j 



with j = 1 . . . 4, and we have omitted the index a for better readability. 

A. 2.1 Partial-wave decomposition of the spin- 3 /2 basis 

In the baryon rest frame, the total spin and relative angular momentum operators are 



S 2 : 
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o^ v (8 1 + perm.) 
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+ 2 L(p) • L (?) , 



(A.22) 



with 



T 2 

r 2 

L (p) ■ L (q) 



2p-V p +p i (p-V p )Vj,-p 2 A p , 
2q-V q + (f(q-V q )V q -q 2 A q , 
P i (q-V p )V q -(p-q)(V p -V q ) , 



(A.23) 



where p and q are the spatial parts of px and qt, respectively. 

It is useful to realize that the basis elements containing {S, V} and {P, A}, which differ 
by a 7 5 (8 7 5 , and those with a different value for the index a = ±, do not mix under the 
action of S 2 or L 2 and then can be analyzed independently. On the other hand, from 
Equations (A.22) and (A.23) one can infer that the set (A.21) can be further subdivided 
into four subsets which (due to their different momentum dependence) do not mix under 
the action of S 2 or L 2 : 
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°11> D 13> °31> 14' D 41> D 32> v 14' v 13 ' 

eg CP qP cP cl cl \tP \tQ 
D 12> D 14' °41' °32' D 13> D 42> v 14' v 13 ' 

c9 cP qP q1 q<? ci \rP \t1 

13 , o u , d 33 , l3 12 , o 34 , o 43 , v n , v 12 , 

C9 Cl Cl CP CP CP \Tl \fP 

14 , a u , o 44 , l3 12 , d 34 , d 43 , v n , v 12 , 



(A.24) 



where the left column indicates symbolically the different momentum dependence of the 
basis elements, in powers of pr and qt, denoted as p and q, respectively. This allows 
to simplify the partial-wave decomposition by looking for S 2 and L 2 eigenfunctions only 
within the above subsets. 

The operator S* 2 is independent of the momentum content of the basis elements. There- 
fore it is sufficient to find the eigenstates at fixed values of pr and % . As we mentioned 
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Table A.l: Orthonormal Dirac basis rf- with s = 3 /2 and for fc = 1. Wc 
omit the Dirac and Lorcntz indices as well as u for better readability. 
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Table A. 2: Orthonormal Dirac basis rf- with s = V 2 and for k = 1. Wc 
omit the Dirac and Lorentz indices as well as a for better readability. 



at the beginning, the problem can be easily implemented and solved using a symbolic 
programming language. 

It is instructive to study how the decomposition for the orbital angular momentum 
is performed using some simple examples (again, we use Mathematica™ for the full 
calculation). The I = elements can be found immediately; they are the momentum 
independent elements in (A. 21), i.e., [Sf-J " and [PfJ' 7 . For the remaining basis elements, 
let us note that they can be written as contractions of 

p a , q a , p a q 13 , pV, ?V, pW, ••• (A.25) 

with appropriate Dirac-Lorentz momentum-independent tensors. For the £ = 1 elements 
it is enough to consider the first three elements in the list above. Applying the orbital 
angular momentum operator one gets 

L 2 p a = 2p a , 

L 2 q a = 2q a , (A.26) 
LW =4p a qP + 2q a pP , 

and from here it is clear that the I? eigenfunctions will come from the combinations 

p a , q a , p a q 13 - q a p 13 , (A.27) 



again contracted with the corresponding Dirac-Lorentz structures. For other £ values the 
calculation proceeds along the same lines but the details are more involved. 
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The above analysis focused on the subset {S,V}, which we will denote by A; = 1. 
Similar results hold for the set {P,A}, denoted by k=2. The final result of the partial- 
wave decomposition is given in tables A.l and A. 2 for k = 1. The case k = 2 is obtained 
from the previous elements by exchanging S — > P, V — > A and adding an extra minus sign 
to the elements ■. This basis fulfills the following orthonormality relation 



l Tr ( =cr,k a',k' 



(A.28) 



/3a8j \ J J a/3j8 

where the conjugation of the basis elements is defined as 

r^ lS (p,q,P) = ~C ah C ld [ T ^ bcd (-p,- q ,-P)] T C^Cj s . (A.29) 

A. 3 Color traces 

To calculate the color factors of the quark DSE and the Faddeev equation we only need 
the following identities for the Gell-Mann color matrices t°"^ B 

N 2 - 1 

= —f^Aii (A.30) 

m,C 

Trt m = (A.31) 

where lower-case and capital indices run over the fundamental and the adjoint represen- 
tations, respectively, of the color SU(N C ) group. 

t 1 

8 to A' 3 A 

i s ^ 

„ e A'B'C x B' B * s ABC 

B t m 

c c 




V AC W Z DB 



Figure A.l: Diagrams to calculate the color factors in the quark DSE (left) 
and in the Faddeev equation (right). 



For the quark propagator DSE we have (see Figure A.l 

r 

" 2N, 



N 2 - 1 

t l AC 5 lm 8cDt 7 B B = ^ t AC t CB = L &AB (A.32) 



l,m,C,D m,C 

and therefore the color factor is (N 2 — 1)/(2N C ), which equals 4/3 for N c = 3. 

For the Faddeev equation we consider the color factor is the same for all three diagrams. 
As shown in Section A. 2, the color part of the Faddeev amplitudes is 

(A.33) 
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Therefore, as depicted in Figure A.l, we have 



E tA'B'C" (,l dm.rn e \ 
JtT- \tA'A d t B'B b C'C ) 



£ABC 

y/6 V™" """"^J^/f 

all indices 



^ ~ CA'B'CCABC (*A'A*B'b) = X/ fi i^A'A^B'B ^A'bSb'A) {^A'A^B'b) 



6 ^ g 

all indices all indices 



1 1 1 JU2 _ i 

El+m j.m j_m + m \ _ \ + m j.m _ c \ y 

g \ l AA l BB - t BA t AB) ~ ~q Z^i T BA z AB ~ ~g 

m,A,B m,A,B c B 

N 2 - 1 



where we used eA'B'C e ABC = ^A'A^B'B — ^A'B^B'A- The color factor thus gives —2/3 for 
iV c = 3. 



A. 4 Numerical details 



A. 4.1 Resolution of the Faddeev equation 

The numerical techniques used in this work are an extension of those explained in [33] 
to the case of the nucleon. We summarize the main ideas in this appendix and extend 
them to the case of the A-baryon. Using the symmetries of the Faddeev amplitudes, and 
requiring certain symmetry properties for the interaction kernel, we will be able to relate 
the three terms in the Faddeev equation (2.26) and, therefore, simplify its resolution. 

The full Faddeev amplitude is the product of color, flavor and spin parts. Since it 
describes baryons in terms of their valence quarks, it must be antisymmetric under the 
exchange of any two of the three quarks, as required by the Pauli principle. Furthermore, 
a baryon must be a color singlet and hence, as explained in Section A. 2, the color part of 
the Faddeev amplitude is cabc an d therefore is always antisymmetric. Thus, the product 
of flavor and spin parts must be symmetric. As we show below, these symmetry properties 
allow to relate the three terms in the covariant Faddeev equation (see Figure 2.6). 

Although in Equation (2.26) (and throughout this work) only the Rainbow-Ladder 
interaction kernel is considered, for the derivation in this section we relax this constraint 
a bit and require the 2-body interaction kernel to be only flavor- and color-independent. 
In this way the color and flavor parts of the Faddeev amplitude factor out and we are left 
with an equation for the spin part 
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^ap~fl(P, q, P) = 



d A k 

(2tt)' 



d A k 



d A k 



Kppi rn 'fa,kz]k) 5 aa ir Spipnfa) Syyifa) 

x ^a"P"j"i(p {1 \q {1 \P) ] + 
K aa > ni ifa, ki; k) Sppn S a i a "(ki) Syyifa) 

x * QT y'i(P (2) '9 (2) . p ) ] + 

Kaa',PP'(kl,k2;k) (5 7 y/ S a ' a "(ki) S p pi fa) 

x Va»p"j"i(p {3 \q {3) ,P) 1 , (A.35) 



where k here only means the exchanged momentum between the two interacting quarks. 

The strategy to relate all terms in this equation is to perform permutations on the 
quark momenta {pi,P2,P3} and quark indices {a, ,$,7} of the Faddeev amplitudes, such 
that the three integrals (which in the following we call II, 12 and 13) look formally the 
same. In particular, for reasons we explain at the end of the section, we want all of them 
to look formally like 13. For the first integral the permutation we are interested in is 
(123) -»• (231) and for the second is (123) -> (312). 

Using now {ABC} as generic quark indices for Dirac, flavor and color indices (e.g. 
A — > {a,a,r} where r would be a color index) and using for clarity {pi,P2,Ps} as the 
arguments of the Faddeev amplitudes instead of the (equivalent) set {p, q, P}, it reads 



Tabcd(puP2,P3)= \^<p 1 i(Pi,P2,P3)®FZ bcd j®^l. (A.36) 

Here the index p denotes the representation of the isospin group to which the baryon 
belongs and F is the flavor part of the amplitude. For example, in the case of baryons 
in the octet SU{3) representation (spin-!/2 baryons like the nucleon), they can belong to 
a mixed-symmetric or a mixed-antisymmetric representation, and the physical state is a 
quantum superposition of both. As explained above, the symmetry properties of the spin 
part can be deduced from the fact that the Faddeev amplitude must be antisymmetric in 
the first three indices. The color part is automatically antisymmetric and this forces the 
product of spin and flavor part to be symmetric. The symmetry properties of the flavor 
parts are known, since they are obtained from the quark model for baryons, and from 
them the transformation properties of the spin part can be inferred. 

The application of this ideas to the nucleon and the Delta have been described in [33] 
and [34], respectively. The nucleon is a member of the octet representation of the SU(3) 
flavor group. There are two such octets, corresponding to mixed-symmetric and mixed- 
antisymmetric flavor states (see e.g. [131]). Under permutations of the quark indices, 
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these two representations mix 



r abcd — IVI 1 r bcad ' 

K bcd = K^Lm , (A.37) 



with p = 1, 2 and 





(A.38) 



and therefore the spin parts must transform in the same way 

*W(Pi,P2,Ps) = Mf '^ x (p 3 ,Pi,P2) , (A.39) 

to get a singlet under permutations for the product of spin and flavor parts. The Delta 
and Omega baryons belong to the decuplet flavor representation, which is a symmetric 
one. In this case, p = 1, Mi 2 = 1 and the spin parts are therefore simply symmetric 
under permutations 

*a/3 7 x(Pl,P2,P3) = ^l3jal(P2,P3,Pl) , 

*a/3 7 x(Pl,P2,P3) = * 7 a0z(P3,Pl,P2) • (A.40) 

In both cases, we show now that the permutation of the quark momenta is equivalent 
to evaluate the amplitude with non-permuted quark momenta (i.e. 13) at different relative 
momenta {p f , q'} and {p", q"}, respectively. First, let us recall the definitions of the relative 
and total momenta 

P 1 - C 

P =(l-C)P3-C(Pl+P2)> Pi =-Q-^ + ^^P, 

P2-P1 P 1 - C D (A 41) 

-P = Pi + P2 + , P3 =P + (P , 

as well as the definitions for the internal relative momenta 



pW = p + k , p^> =p — k , = p , 

qW=q-k/2, qW=q-k/2, q^ = q + k , 



(A.42) 



with fcj = pi — k and fcj = pi + k the internal quark momenta. Now, let us show that 
the quark and internal relative momenta in the terms II and 12 of the Faddeev equation 
can be written as those in 73 by introducing new external relative momenta {p',q'} and 
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{p",q"}, 



13 : < 



+ 



k 2 = q~l + ^ + k 
[ P3 = P + f 
(Pi = ~Q ~ \ + 



II : < 



(Ps =P' + 



(k 2 = q - \ + f - k) = (k[ = -q' - £ + f - k) 
[ (fc 3 = p +f + fc) = (# 2 + f+fc) 



i p i q . 3p 



(A.43) 



12 : < 



(%! = -q - I + f + k) = (k 2 ' = q" - ^ + f + fc) 
[ (k 3 = p + ^ - k) = {k'l = -q" — ^- + f — &) 



P =9-2 > 



9 2 4 



(A.44) 



This identification is possible if one has (1 — Q/2 = £ in (A. 41). To this end, we need to 
set C = 1/3. 

Making these substitutions in (A. 35), we get 



+ 



^ 2 ' (27T)4 



d 4 k 

J2ny 



Kpp'^'ik'^k^k) 5 aa » Spipi^k'-t) Syy>(k' 2 ) 
x* p ;,, ra ,, x (p^,q^,P) 

K aa ',jj'(ki,k 2 ,k) Sppi Syj"(ki) S a > a "(k 2 ) 
x $> p ' a"W P) 

W j"a"/3"X\P 'V i r ) 

K aa >,/3i3'(ki,k 2 ,k) (5 77 » S a ' a "(ki) Sppn(k 2 ) 



+ 



with the internal relative momenta defined, as in (2.30), 
pi 1 ) = p-\- k, p( 2 ) = p — k, 



p( 3 ) = p : 



qW = q - k/2, qW = q - k/2, q^ = q + k. 
Renaming dummy indices it becomes clear that, if the interaction kernel is such that 



(A.46) 



K a a',-yY {k\ ,k 2 ,k) — Kjji aa i (k 1 ,k 2 ,k) 



(A.47) 
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then the three integrals are formally the same and, denoting 13 by ^^^(p, q, P), we have 

so that in practice we need to calculate only Vl/( 3 ) and by further evaluating it at the points 
{p', q'} and {p", q"} we obtain the full Faddeev amplitude. 

The requirement (A. 47) is obviously satisfied by the Rainbow-Ladder kernel (2.21) 



K 



a eff(k 2 
k 2 



V 

aa' '77' 



(A.49) 



Nevertheless, the condition (A. 47) does not seem to be a very stringent one and it may 
be fulfilled by many 'beyond Rainbow-Ladder' kernels. 

In practical calculations one solves for the scalar coefficients in the expansion of the 
Faddeev amplitudes 

9. P) = & m (p 2 , q\ zo, Z!,z 2 ) T^ip, q, P) , (A.50) 
with zo = pr ■ qr, Z\ =p-P and z 2 = q- P. Then, the Faddeev equation for the coefficients 

/ W,«)(p2 > g 2 > ZOj Zu Z2) =/ (p),«).(3)( p 2 j g 2 ) 2()) Zlj Z2 ) + 

Mr'^/^>0),(3) (p /2 ) ^2 )4)4)4) + 



with 



H[ 3 



TT13 

n 2 



^ a iy(P,q,P)T^ aX (p',q',P) 
f l p aXj (p,q,Py iapi (p",q",P) 



(A.51) 



(A.52) 



and 

/ W ' W ' (3 V>«Wi,*a 



d A k 



(2k) 



Tr 



T jfal7 (P> P ) K ota'W (h ,k 2 ,k) 5-, 
X W(fcl) ty^fe) T$M^\<P\P) ] f^ m3) ((p^) 2 , (9 (3) ) 2 , 4 3) , 4 3) ) • 

(A.53) 



Let us show why it is convenient to represent all diagrams in terms of 13. According to 
(A. 46), in this case the internal and external relative momentum p is the same (p^ = p), 
and can be chosen as in (A. 3). On the other hand, if we change the integration variable* 
from the gluon momentum k to the internal relative momentum (the Jacobian of this 



"This is, in principle, only possible if we use a translation-invariant regularization of the integrals. 
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transformation is, according to (A. 46), unity), then the transverse momenta which appear 
in the definition of the Faddeev basis (see Section A. 2) simplify to 



PT 



(3) 
P T 



/o\ 



1 

v o y 



qt 



l 



V o y 



/ 



(3) 



V 



sin ip \ 
cos ip 


) 



(A.54) 



with <f an angular integration variable. Now, defining the Faddeev wave-function as 

$a/3 7 x(P, q, P) = S aa , (pi) ,%/ (P2) ^a'P'-y'liP, 9, -P) (A. 55) 

it can be expanded in terms of the same Faddeev basis. Then, the interaction kernel 
matrix becomes 



K 



(3) 



(*0 



^(P, 9, P) [iLnh'hi') <^ YI (p (3 \q {3 \P) 



(A.56) 



and depends only on the angular variable p>. In the numerical resolution of the Faddeev 
equation, this kernel is treated as a matrix, which is now small enough to be calculated 
in advance, stored in memory and reused during the iteration process. This reduces 
considerably the computation time. 

Finally, let us briefly describe how the Faddeev equation (A. 51) is solved numerically. 
The first step is to modify (A. 53) by a multiplicative factor A, thus transforming it into an 
eigenvalue equation. The physical solution corresponds to the case A = 1. The simplest 
method to solve an eigenvalue problem is by iteration. A test baryon-mass as well as a 
starting function for the coefficients /(^'W in the integral (A. 53) must be chosen for the 
first iteration. The resulting functions are used as an input for the next iteration and so 
on. If the converged eigenvalue is not 1, then one chooses a different value for the test 
mass and repeats the procedure. 

To optimize the numerics, we expand the angular dependence of the coefficients in 
Chebyshev polynomials 



f^(p 2 ,q 2 ,z ,z 1 ,z 2 ) = fj±®tfrf )Cl fa )Cn{zi)cM 
with C n (z) = i n U n (z) and U n (z) the Chebyshev polynomials of the second kind 

7rk 



Un (z) = 2 n f[(z- cos , (n>0) 
fc=i ^ n ' 



(A.57) 



(A.58) 



The main advantage of this expansion is that the angular dependence of the Faddeev 
amplitudes is typically weak, and therefore a very small number of Chebyshev moments 
finm^^P 2 ^ 2 ) are nee ded (for an illustration of a typical case, see Figure A. 2). 

Although the resulting baryon masses are not very sensitive to the precision used for the 
numerical integration in the Faddeev equation, for the calculation of form factors one needs 



A.4. NUMERICAL DETAILS 71 




-0.2 1 1 ' 1 

0.5 1 1.5 

p(GeV) 

Figure A. 2: Chebyshev moments of the z\ variable as a function of p, for 
q = 0, zq = and z-i = 0. 



a very precise determination of the Faddeev amplitudes. For the calculation in this work 
we used 30 Gauss-Legendre quadrature points for {p 2 ,q 2 }, 8 Gauss-Chebyshev points 
for {z, zo, z±, Z2} and 8 Gauss-Legendre points for {y,<f} (using the notation of Section 
A.l). We use four Chebyshev moments in the expansion (A. 57). For the evaluation of the 

amplitudes fi^Jf' \p 2 ,q 2 ) at {p 2 , g 2 }-values different from the Gaussian quadrature points, 
we use cubic-spline interpolation. 



A. 4. 2 Form factor calculation 



Using the same transformations we introduced in previous section, the three terms in 
the equation for the electromagnetic current in Rainbow-Ladder truncation (3.19) can be 
expressed in terms of one of them. In this section, we limit the discussion to the case of 
spin-3/2 baryons (the spin-!/2 case has been treated in [33]). 



72 



APPENDIX A. 



The equation for the current is 

J T>T= Gl / / ^/3'a'Xv(P/ 1} >4 1} ' P /) [( 5 (Pl) r ' 1 (Pl>Q) 5 (M)) Q/Q ^ /3 (p2)Sy 7 ( K 

+ Q 2 y j^^TY(p { f 2} ,qf } ,Pf) S a/a ( Pl )( y S(p f 2 )T^p 2 ,Q)S(pl))^^S Y7 (p 3 ) x 

(^ 7 x(pf } ,,f } ,PO-*g 7 x(^ } ^f } ,^)) 
+Q 3 | J^(3'a'x^(pf } ,q { f 3} ,P f ) S ala (p 1 )S (S ^(p 2 )^S(pi)T^ P3 ,Q)S(pl))^ x 

(*a^(p{ 3} , rf 3> , Pi) " * S } 7 X(PP } > ^ 3} > P)) - 

(A.59) 

where we now wrote explicitly the charge of the quark i, Qi, that is part of the quark- 
photon vertex. Let us apply the transformations (A.40-A.44) to this equation. For exam- 
ple, the first term becomes 



Qi 



%^iw( P f } ,q' f {3} ,Pf) (S(p'l)r»(p' 3 ,Q)S(pl)) / SppbtiSfM) 



pJq 



(*W(p? 3} , ,?\ Pi) - *&r(p? 3} , ^ Pi)) (A.60) 
where we used the definitions (A. 43) and (3.22), with £ = 1/3, and we also introduced 



/ P'P pP 
P ' = - q ~ 2 + 3 =q ~2 + 3 =P2 ' 

/ / P' P P 

P2 = Q ~^ + J = p+ ^ =P3 > 

i i P P P 

P 3 =P+^ = - q - 2 + j=Pi. 



(A.61) 



Performing the change of variables {p, q} — > {p', q'} in the integral (which is a transfor- 
mation of Jacobian unity), this term is exactly equal to the third term in (A.59). In a 
similar way, the second term in (A.59) can be written formally like the third one. There- 
fore, introducing the total baryon charge Qb = Qi + Q2 + Q3, the calculation of the 
electromagnetic current simplifies to 



4, T = Qb 



pJq 



ypwi>y(pf\qf\Pf) 



ScMSpM (s(p()r»( P3 , Q)s(pi)) 

(<W(P? } , ^ {3} , Pi) - * { X x {p?\q?\ Pi)) , 

(A.62) 

and it is clear that for a neutral spin-3/2 baryon, all electromagnetic form factors will vanish 
identically. This as a drawback of assuming isospin symmetry in a covariant Bethe-Salpeter 
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approach. Finally, this simplification is valid only in Rainbow-Ladder or, more generally, 
whenever the interaction kernel does not couple to the external photon and fulfills (A. 47). 

As already mentioned in previous section, a reliable calculation of (A. 62) requires a 
very precise determination of A naive approach to solve (A. 62) with the number of 
integration points used in this work, implies an evaluation of (128 x4x4x4)~6 - 10 7 
functions (each of the coefficients in the expansion (A. 57)) at (30 x 8 x 8 x 8) 2 ~ 8 • 10 8 
quadrature points. It is clear that one needs to arrange the integration in a way such that 
the number of evaluations is reduced. In this respect, if we use the Faddeev wave functions 
(A. 55) we can write Equation (A. 62) as 




x 



(w( P { 3, , g { s, ,p*)-*a I (p{ s} ,(rf s} ,po) 



(A.63) 



and since p% and 



j>3 are independent of q, we can distribute the integration 




x |*/^y(pf,ffj 3 ^)(^^ • (A-64) 



Moreover, it is easy to check that the momenta qJ/J, that appear in the arguments of the 

Faddeev basis, and the angles {z l J^ , , z^} are independent of q 2 , which simplifies the 
numerical integration. 
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